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INTRODUCTION 


In  parallel  with  a  continuing  active  use  of  the  AFGL  Rocket  Trajectory 

System,  improvements  and  extended  features  have  been  incorporated  in  nearly 

all  phases  of  this  system.  This  combined  report  and  user's  guide  documents 

current  procedures  and  provides  the  analytical  background  underlying  develop- 

1  ?  1 

ments  implemented  since  the  installation  of  the  original  system.  ’ 

The  user  is  referred  to  the  above  referenced  reports  for  details  of  the 
scope,  analysis  and  implementation  of  these  programs.  Major  revisions 
resulting  in  the  present  trajectory  system  include: 

1)  A  combined  editor  and  preprocessor  program  DRIVEA,  which 
replaces  the  original  DRIVE1-DRIVE2. 

2)  A  comprehensive  filtering  program  DRIVEB,  which  permits 
launch  optimization  estimation  and  segmented  integration¬ 
filtering-reproducing  modes. 

3)  Substantial  augmentation  of  TAPE4,  the  DRIVEB  output  data 
file,  to  include  residual,  launcher  referenced,  ballistic, 
and  local  velocity  information. 

4)  A  simplified  multi -radar  solving  program  DRIVEC,  which  is 
based  on  more  reliable  estimation  of  the  relative  error 
covariance  matrices. 

5)  Streamlining  of  the  processing  system  by  elimination  of 
buffering  and  CPRS,  elimination  of  the  DRI VE1-DR IVE2  pass, 
reduction  and  flexible  interpretation  of  card  data  input, 
tape  processing  aids  in  DRIVEA,  standardization  of  common 
routines,  and  the  use  of  CDC  utilities. 

This  report  is  written  so  as  to  be  maintained  in  ring-binder  form  for  easy 
revision  and  updating.  Section  2  presents  the  input  card  data  formats  and 
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the  tape  Input  and  output  formats.  This  information  is  also  commented  into 
each  of  the  programs.  A  brief  description  of  the  various  modes  available 
in  0RIVE8  is  also  presented. 

Section  3  outlines  the  organization  of  each  program.  The  supporting  material, 
which  either  analyzes  new  developments  or  attempts  to  substantiate  original 
procedures,  is  presented  in  the  Appendices. 
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2. 


USAGE  GUIDE 


The  AFGL  Rocket  Trajectory  System  is  designed  to  accept  raw  radar  data  from 
any  rocket  range  on  magnetic  tapes  in  the  format  shown  in  Figure  1.  It 
then  generates  an  intermediate  edited  file  and  associated  print-out  using 
program  DRIVEA.  Finally,  it  provides  a  filtered  best-estimate  trajectory 
with  accompanying  reports,  plots  and  summary  file  using  program  DRIVEB  in 
one  of  a  variety  of  optional  filtering  modes.  In  the  event  that  multi-radar 
tracking  data  is  available,  a  statistically  weighted  final  trajectory  with 
reports,  plots,  and  summary  file  can  be  developed  using  DRIVEC.  The  user 
outputs  from  DRIVEB  and  DRIVEC  are  identical,  and  the  summary  file  (TAPE4)  can 
be  used  at  any  time  to  recreate  the  reports  and  plots  using  program  REGEN. 
These  four  programs  comprise  the  AFGL  Rocket  Trajectory  System.  In  this 
section  operational  modes  and  input  and  output  file  formats  are  outlined 

for  each  program.  The  output  files  with  their  headers  also  serve  as  the  input 
for  the  subsequent  programs. 

2.1  DRIVEA 


Figure  2  presents  the  card  input  format.  Details  of  calibration  and 
refraction  correction  are  unchanged  from  Reference  .  A  number  of  improve¬ 
ments  have  been  made  over  the  original  DRIVE1-DRIVE2  system.  T0L(2)  and 
T0L(3)  control  data  tolerance  and  modification  (25%)  limits  are  described 
in  Appendix  A.  NSKIP  controls  initial  positioning  of  the  raw  data  file, 
and  NOK  indicates  which  samples  in  the  next  record  are  to  be  used  as  the 
starting  base.  TLO  and  THI  control  the  first  and  last  data  points  to  be 
entered  on  the  edited  output  file  (Figure  3)*  the  actual  times  are  recorded 
on  the  header  as  TSTART  and  TSTOP,  and  are  later  used  automatically  by 
DRIVEB. 


Spanning  of  Data  Gaps 

An  additional  feature  available  in  DRIVEA  allows  spanning  of  missing  or 
poor  data  by  inputting  VAR,  TI,  and  TF  for  smooth  polynomial  substitution 
of  range,  azimuth,  and/or  elevation  samples.  VAR  R,  A,  or  E  identifies 
the  variable  to  be  replaced  over  the  time  span  TI,  TF.  Up  to  two  cases  are 
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TAPE1  -  RAWDA.T 


This  Cape  is  a  standard  raw  data  format  and  consists  of  a 
20-word  header  followed  by  a  series  of  96-word  data  blocks.  Each  data 
block  contains  6  complete  data  samples  consisting  of  16  words  describing 
the  sample.  The  editor  does  not  use  all  of  these  words,  and  only  those 
shown  below  are  required. 


HEADER  RECORD 


WORD 

MODE 

VARIABLE 

2 

Floating 

Radar  latitude  in  radians. 

3 

Floating 

Radar  longitude  in  radians. 

A 

Floating 

Radar  height  in  feet. 

5 

Floating 

Direction  of  flight  line  in  radians. 

8 

Integer 

Launch  day. 

9 

Integer 

Launch  month. 

10 

Integer 

Launch  year. 

DATA  RECORD 

WORD 

MODE 

VARIABLE 

1 

Floating 

Time  U.T.  in  milliseconds. 

6 

Floating 

Slant  range  in  yards. 

9 

Floating 

Azimuth  in  degrees. 

10 

Floating 

Elevation  in  degrees. 

| 


Figure  1  Format  of  Raw  Range  Radar  Data--Tape  1--Input  for  DRIVEA 


PROGRAM  CRIVEa<INPUT,eUTPUT,ftAW0ATtTAPEll,TAOEl=»AWCAT 
—  +»-TAPet<^,-T-AF€9»-YiAPFB--,TAPE«r=OtMYFT.TT-,-YA»E-5*TM«ttT) - 

-optv»A — we-Ygccpga  by  ttsieett-.iwe. - t^w-ra* — — 

•••••  *•»«* 


i 

( 

! 


c 


I  1 


THIS  PROGRAM  COMBINES  FUNCTIONS  PERFORMED  BY  TRAJECTORY  ANALYSIS 
PP£)6RAMS~»MTfY«--»N«--9friYe-?-.  — TH?-»YMA^C  FlLTER— CF — DR-iVc-E- H-AR  BEEN 


REPLACED  3Y  A  OATA  MCD IFICATI CN/ I NTERPOLA TI CN  FbCCSDURE  IN  THE  EDITOR 


P9EFR<3CESSOR(,’OIYS2>  PHASE,  THREE  PARAMETERS  (T  r  l  ( 2)  ,Tct  (3 )  ,  NOK> 

-HA**-  YE  CM  ADDED  FOR  THC  POINT  HBCIFICATIO  M/t  NT  CPR-Et  AT-'HON-  PRDO-eOUcP-t - 

BUT  DEFAULT  VALUES  ARE  AVAILABLE.  CRIVEA  RESULT*  IN  FASTER 


0RIVE2  FILTER  ( I  •  £•  ff  IL  TER  RINGING). 


-eaae  hypuf  -rewirementRi 


NO.  VAPIAELS 

-  1 - 

PLIGHT 


DESCRIPTION 

FLIGHT  IDENTIFICATION,  20  CHARS. 
-CHARACTERS-,  THESE  ITEMS-  wHrL - 


APPEAR  on  THE  PRINTED  OUTPUT. 
-a-9-  gH-ARACTER"  Y-I-tlE  FDR-  T^ft-THCR 


-4-Af3~ 


RUN  IOENTIFICATICN.  MILL  APPEAR 

- e-H  F  c  I>  TO  O  OUTPUT  -CNfc-*-. - 

I COO E  RANGE  IOENOIFICATICN  code. 

- 1— FflF — EGL-I  N--A.  F-rC-» - 


15 


2  FOR  NALLOPS  IS. 


hfi  ; 


OATA  SAMPLING  TIME  INTERVAL 
-FN- 

TfftX ,  MAXIMUM  EDITING  TIME  IN  SECONDS. 

MEASURED  AS  a  DISPLACEMENT  p=cm 

- ___ - -THg  FIRST  OATA-FOOT  ON  RA-MOAT,- — - - 

“rOL  (  2  J  ,  TOLERANCE  LIMIT-  DEFAULT  VALUE  (3.0)* 


,215 


NO*, 


POINTS  6-NO<  I  5-NOK  ARE  TAKEN  AS  VA LI r 
»«Sg"FCIMTS-  IN  TH*  FIRST  OATA  0« 
NUMBER  OF  INITIAL  SUFFER  RECORDS  TO  °E 


SAMPLE 
FP  EO. 


N**IP 


p<: 


NOTE,  the  FOLLOWING  T  OATA  CAROS  ARE  NOT  REQUIRED  UNLE'S  T“E 
- - <RAHGE-~CCCC~  SFECI'FIES-EHURCHIL'L. — <-fvE» — fCCDE^T) - 

* - NCODEi - Dt/OE-TH  DETERMINE1  -  IF-  A  RANGE  -COR- IB;EF  13. 5 - — - 

RECtICN  IS  REQUIRED. 

- * — f -f-OR  PA  F6E  CORRECF^eHr— 1 RAfrAR— NQ-,  1  >■■ -»?  -RAOAF  HO.E-.~ 

*  0  SETS  RAN"  CCH SECTION  TO  TIRO. 

- i-f  nc-corrcc  tion  is  »;wrw  th-7 - 

REMAINING  FIELDS  NAT  RE  EEROS, 

- - CAL-RBi - R*NG€-CAt-I<?R*MON-  CISTANOE-  IN  YAP-OS, - 

cs,  RANGE  CALIBRATION  READING  IN  YARCS. 

-P*S - *€•*€£*■  PEL  AY  IN  Mir RR3*E-G-6 NO B-« - 

RTeM ,  REACCn  TRANSMITTER  PULSE  WT0th  j k 

- - — - 

PTFM  RADAR  TRANMITFR  oilSE  WIDTH  JN 

- HieReDcceMTOi - - - 


Figure  2a-  Card  Input  Format  for  DRIVEA 


ncooe, 

4Z=, 

A?* - 

i  *  VJ  U  U  L  | 

*LP, 
ELN 

3  KSVQPT 


1**3-, 

ITB, 

Wr- 

TMfS, 


T^EC 


*1  FOR  A2IHUTHAL  COPRECTION.  IS,2F15.3 

■3  SSTE-  iPI-HUTH  GORREOT£GN-TG—  *G3», - - 

A  7TMUTN  °LUNG£  READING  IN  CEGREES. 

A MHUTH-  BEADING  IN-NCRHAL  POSITION - 

IN  0?GeEES. 

.»<  "FOB  -etE-YATION-eefi^COTIONi - IS,  gFi  3,3 - 

*3  SETS  ELEVATION  CORRECTION  to 

ZERO.  - - 

ELEVATION  REAOING  IN  THE  BLUNGE  =0- 

SITICN  -IN— OE GR5SS. - 

ELEVATION  REAOING  IN  THE  NORMAL  FO- 

■9 Wf  0 N  -tN-RS TRCrr, - 

OPTION  KEY.  15 

H- FGR-NE-GEOGENTRIG  CONVERSION. - 

»2  FOR  GEOCENTRIC  CONVERSION. 

C  ALENO  AB-  SRT E4H6WTA  , OAT,  TE A B> - M-I-s,  fr*4-r IV,3F-l-» 

AND  U.t. (HOURS, HINLTES, SECONDS)  OF  LAUNCH. 


-HEIGHT - HE  IGHT  -OP- c A  OAR  A90VC  TltC - F-JrS-.-e - 

SPHE»0I0  IN  “Etess. 

■ftWeHy - RADAR  tONGIT L*1C  IN1  BCGRClO , HI N i  , - SPfrfrrN - 

CLCMIN,  ANO  SECONDS.  NEGATIVE  IF  EAST. 

'BfeC'RrEO  - - - — — - 

RLAOEG,  RADAR  I.ATITUOE  IN  CEG=EE  S ,  H  IN  . ,  3R13.A 

-H.AHIN, - AN"  SECONDS. - NEGATIVE  IF  SOUTH. - 

'LAS  EC 

*3 - FfcEr - T  INE  Ll-H**R -B4-  SSCGN85- U»  t.  pOB - iTHri - 

TMI  GEOCENTRIC  CONVERSION. 

-  DEFAULT  VALUES  (O.OtV - - 

TLOaTlHE  of  FIRST  ran  "ATA  POINT. 

THI*TIHE  0p  LAST  RAH  DATA  POINT  PROCESSES 
- ar  bROGRAH  . - 


IN 

IREF 

REFRACTION  CORRECTION  <£v.  1= 

REFRACTION  needed, I. E.,  the  "'ATA 

15 

FPOGIO 

- T3 — WCW^wr^i — Ll  1  .'Jr - 

three  voph  pROGpan  ioentifie* 

1502  EJnsSC  bbIVA 

V  *4  “  lUDkL  1  t  ini  n  -<  if  .  J  *  J 

SHCCTH  P0LyNG1I A L  ^INSTITUTION  1$  "ESI^FC. 

rtf 

- *e - 

START  ANO  Esc  OF  T IMES  FROH  LAUNCH  ISEC'.) 

TAPE  11  IS  A  SEQUENTIAL  FILE  CF  th£  EDITfo  ^sh  1JTA,  An3  SERVES  r  A 

~HNR  *STNfgN  DRIVB1  ANfr  DRtVS*  PHASED. - TARES -fry  St  *Nfr"H - A-R* - 

STNAOv  fIL£«  HOLDING  the  final  GEOCENTRIC  VELOCITIES,  GECCENTsi: 
COORDtHAVEB  I  AND  EDITCC  RAW  •R-ATA-  RESPECTIVELY-, - - - 


FIGURE  2b.  CARD  INPUT  FOR  DRIVEA 
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* 


TAPE 

10  INCLUDES 

A  HEADER  RECORD. 

woscOfST-svrteet 

-UESCRI-PTI-CN-— 

-  .  ...  . 

1-2 

OATS 

RUN  DATE,  BLANK 

3-4 

— DAT  tJ-4) 

-  FLIGH-T-fOtfmFI CATION- . 

•  —  - 

5 

LOAT 

LAUNCH  CATE  (MO/OA/YEAR) 

6 

—  TLNGH - 

-LA UMG F—TTME  -UT  {  S  50  S> 

. 

7-9 

PROGinCl 

-3PR0GRAM  10  FOP  PLOTTED  OUTPUT 

1C 

-Hf-AO - 

-RADAR- AfcTTTUEHE  (KM) - 

11 

RLOCSG 

PAOAP  LCNCITUOE  (CEG) 

12 

PJLGHIM- 

—  ’RADAR- LCNGItU0E  (MIM)  -  - 

- - 

13 

RLOSEC 

RADAR  LONGITUDE  (SEC) 

14 

P.LAC2C- 

RACAR  LATITUDE  (OEG> 

.  ........ 

15 

*?LA  PIN 

RAOAP.  LATITUOE  (MIN> 

It 

--  RlfrSEC  - 

-  RADAR  LATITUDE  (SEC) 

... 

17 

TSTART 

DATA  STARTING  TIME  UT  (SECS) 

ia — 

—  —  "^ST€F‘  — 

-  FT  N  AL  •  D  AT  A-  TITHE-  IFF— C-SE&Sf - 

. . — 

19 

OT 

SAMPLING  TIME  INTERVAL  (SEC) 

20-  - 

- JO  O  LE - 

— i  FOR— EGLTN~>-  2  FOR- WALLOP^, 

?  FOR  CHURCHILL 

21 

IFEF 

IF  .GT.  0  OATA  HAS  CORRECTED 

FCR  REFRACTION 

FOLLOWED  BY  RANGE,  AZ,  EL  (KM, RAD)  LOGICAL  RECORDS. 


i- 

r 


Figure  3  ' 

FORMAT  OF  DRIVEA  OUTPUT  FILE,  TAPE  10.  ii 
TKTS  IS  A  BINARY  FILE  GENERATED  BY  DRIVEA  j; 
WH'CH  HOLDS  THE  EDITED  RAW  DATA  USED  AS  INPUT  i, 
FUR  DRIVEB.  jj 


-13- 


allowed  for  each.  Good  samples  at  TI-5,  TI,  TF  and  TF+5  seconds  are  used 
to  define  the  spanning  cubic.  If  samples  are  unavailable  at  TI-5  or  at 
TF+5,  all  times  are  readjusted  so  that  at  least  2-second  separations  are 
maintained.  Special  care  is  required  in  selecting  each  TI,  TF  pair,  as 
described  in  Appendix  A. 3. 

2.2  DR1VEB 


Before  the  most  satisfactory  best-estimate  trajectory  can  be  determined, 
the  multi-option  filtering  program,  DRIVEB,  is  usually  run  on  the  edited 
radar  data  file,  TAPE  10,  for  a  variety  of  choices  of  filtering  time  limits 
and  radar  parameters.  Figure  4  presents  the  card  input  format. 

Radar  Parameters 

Radar  parameters  are  generally  prescribed  for  any  given  radar,  but  the 
investigator's  knowledge  of  the  evident  quality  of  the  raw  radar  data  may 
be  imparted  through  adjusted  values  (see  Appendix  G)  to  more  appropriately 
influence  the  statistical  filtering  equations. 

Ballistic  Parameters 

Rocket  drag,  thrust  factor  and  delay  time  may  be  input  to  override  nominal 
values  which  are  built  into  rocket  models  in  the  THRUST  routine.  Rocket 
models  presently  included  are  shown  in  Table  1;  the  key  is  the  identification 
number  contained  in  DAT ( 5 ) .  Interpretation  of  filter  time  limits  TB...TX, 
KEYOPT,  and  the  geocentric  position  and  velocity  Pl(l-6),  are  described 
with  the  aid  of  Figure  5. 


Run  Options 

Except  in  the  pure  integration  mode  (KEY0PT=0),  TSTART  is  redefined  by  the 
TAPE10  header.  The  type  of  DRIVEB  run  determined  by  KEYOPT  refers  to  the  main 
filtering  segment  TC  to  TE.  The  geocentric  position  and  velocity  values  refer 
to  time  TC  and  may  initially  be  obtained  as  the  rough  estimates  from  the 
DRIVEA  print-out.  The  exception  to  this  is  for  KEY0PT=1  or  -I  when  launch 
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PROGRAM  OR! VE8< INPUT,  QUTPUT»TAPGi(l*TAPE4,TAPE7 
1,  r4PE6*0UTPUT,rAPE5*IM»Ur> 


DRIVER  OEVELOPEO  BY  LOGICON.INC.  Y-P5ION  SEP.,  J')T4 

**»*♦ .  ....  •**»• 

I/O  Dcc I NITIONS . 

T4»E10  RAW  DATA  INPUT  TAPE,  INCLUDING  HEADER. 

this  tape  is  generated  dy  dpivea. 

T APE 4-  PRINTED  ANO  PLOTTED  OUTPUT  UAT A  CONTAINS 

SELECTED  POINTS  IF  TOPT*0,  ALL  POINTS  TF  r0PT=1 
TAPF7-.  FILTERED- AND- RAH.  DATA  OUT=JT-HASS  STORAGE  ACCESS 
ORGANIZATION  FOR  FILE  MANIPULATION 
TAPES-  SYSTEM  PRINTER. 

T 1°F5  SVSTEH  CARO  RtAOER. 

ORIVSR  FILTER  CARD  INPUT  DATA. 

F  1C  21  5* 

CAPO  ITEM-  -  -  DESCRIPTION.  FORMAT 

DATE  ANO  TITLE  INFOPMATI ON.  AA1C 

1  TTTLcl  ANY  SO-  CHARACTER  TITLE 

2  TITLE2  FIRST  20  ANO  LAST  20  CHARACTERS  ART 

IGNORED-  THE  REMAINING  AO 
CHARACTERS  MAY  3E  USED  rOP 
FURTHER  I 1ENTIFICATI ON . 

3  DAT (3-6)  LAUNCH  DATE  ANO  ROC<ET  NUMBER  5  USE  UP  T" 

»»  CHARACTERS  EACH. 

T NO  DIGITS  IN  COLUMNS  21 -2°  FOLLOWED  BY 
A  PER IOO  DETERMINE  rocket  type; 

IF  THE  PER  1 00  IS  ABSENT,  AN  ARBIYR  fRv 
UNPOHEREO  ROCHFT-IS  ASSJMEO. 

4  ooOGIO  PROGRAM  ID  FOR  PLOTTED  DUTPUT.  3P 

CHARACTERS-  FIRST  5  ARE  THE  ID  FRO  * 

THE  JOB  CARO,  NEYT  4  ARE  ThE  CHARG c 
NU-  REMAINING  CAN  PE  JSE 0  FOR  FURTt-EP 
IOENTTFICA  TTON. 

-  IE-  CSXHGi-15  02-3CNO  N  BE  IN.FEB.  2fl»  72T 
RAOAR  POSITION  OATA.  3^10.4 

5  MRAOM  HEIGHT  OF  THE  RADAR  A  3 ON E  THE 

SPHEROID  IN  HETERS. 

6-  RL00E5—  LONGITUDE  IN  DEGREES,  MINUTES 

RLOMIN,  ANO  SECONDS.  IF  THE  .DNCITUDE 

...  RLOSEC -  I*.  EAST  ENTER.  AT  LEAST  DNELOF 

THE  QUANT  I TIES  AS  A  NEGATIVE 
NUM9ER- 

7  PLAJEG,  LATITUDE  IN  DEGREES,  1TNUTES 

RLAMIN,  ANO  SECONDS.  IF  THE  -ATI TUOE 
RLASEC  IS  SOUTH  ENTER  AT  LEAST  ONE 

_  .  D£  THE-  QUA.NTITITES  AS  A  NEGATIVE 

NUMBER. 

RAOAR  PARANETERS.  cc10.4 

A  RPVU),  PEAK  TRANSMITTED  POWER  IN  KILO- 

HATTS- 

R°V (2)  ,  ANTENNA  GAIN  IN  3D.  (OOeS  NOT 
-  ...  INCLUOE-  ANTENNA  EFFICIENCY)  _ 

RPV (3) ,  TRANSMITTED  WAVFLEN5TH  IN 
METERS. 

RPV ( 41 ,  RECEIVER  NOISE  FIGURE  IN  OB. 

RPV  ( 5  >.  SKIN  TRACK  LOSS  FACTOR  IN  OB. 


Figure  4a.  Card  Input  Format  for  DRIVEB 
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9  RP V(6).  BEACON  TRACK  LOSS  FACTOR  IN  03. 

Rpvm,  angle  velocity  la-;  cjefficient 

IN-  MlLS/HtLS/SEC. 

RPV(8>,  RANGE  VELOCITY  LAG  COEFFICIENT 

-  - -  _  -*«S*lC*OS/SEC.-  -  -  .  _ 

BPV(9>,  ScACON  POWER  IN  WATTS. 

RPV(1C>-  BEACON  ANTENNA  GAIN  IN  03. 

f*  RPV(11>,  ANTENNA  8:  ANWIfTTM  IN  DEGREES . 

?PV(U).  ANGLE  TRACKING  »£RVD  4 ANOMIOTH 
IN  HZ. 

-RPVtlSJ..  TARGET  -CROSS  SECTION  IN  03. 

R®V  (14)  i  9ASE  ANGULAR  VAP.IAnSE  IN  WIIS*»2. 

RPV  ( 151  BASE  RANGE  VARIANCE  IN  YAROS»»Z. 
tl  R  PV ( 161 »  RANGE  TRACKING  SERVO  9 ANGM IOTH  IN  ►  ». 

RPVC17L.  PULSE  LENi-TH  IN.  HICROSEC ONOS. 

7PV  (18  (  ,  I.F.  3AN0WIDTH  IN  MHZ. 

R  PV  <  19)  .  TINE  OF  LOSS-  OF  BEACON  TRACK  SIGNAL 

IN  SECQNOA  U.T. 

R  PV  (20)  CONICAL  SCANIGT. 0.81  OR  NONOPULSEC  >0.0) 

ROCKET  ’AYLOAO  AND  PRAG  DATA.  5F10.4 

1Z  * AYLD.  PAYLOAD  NASS  IN  KILOGRAMS. 

RAREAM,  ORAG  CROSS  SECTION  I A  N£TtRS*»Z. 

■If  ZERO,- JSE  VALUE  F'ON  ROCKET  TABLE. 

COINP,  ORAG  COEF*  ICIENT  AT  RACK  VELOCITY  *. 

IF  ZERO.  USE  VALUE  cROM  ROCKET  T A® LE* 

THFACT,  THRUST  FACTOR  ( P£  FAJLT ■N0MINAL*1. ) 

OFACT  OELAJ  FACTOR.  S»»*T  0s  THRUST  IN  Sf"P. 

AFTER  LAUNCH.  (SHOULD  NOT  BE  NEGAT  1VE1 
LAUNCH-  OATS  AND-II-MF.  —  -  -  Z  ( TZ  ,1X)  .14  .3F10.4 

13  LMO,  MONTH, 

LOA.  DAY, 

LYR,  YEAR. 

JHSS.  UNIVERSAL  TINE  AT  LAUNCH  IN 

JMIN,  HOURS,  MINUTES  ANO  SECONOS. 

nsec--  -  —  -  - . 

DATA  STARTING  TIME  ANO  SAMPLING  n  "E  INTERVAL.  4^10.4 

14  THRS,  TIHE.-JT0,  OF  THE  FIRST  DATA  POINT 

TMTN,  ON  TAPE10  IN  HOURS,  HHjrES, 

TSEC.  AND  SECONDS.  DT  IS  HE  TIME  INCREMENT 

OT  BETWEEN  SAMPLES. 

NOTE.  -TIME-  L  OT-- FROM- CARO- 14--ARE  USES'-  ONLY  FOR  KEYO1»T=0, 

OTHERWISE  THEY  ARE  08TAINEO  FROM  THE  TAPE18  VEAPPR. 

I  NT  CGRA  TION.  FILTERING  ANO-  REPRODUCING  TIME  LIMITS.  6F10.4 

15  TP,  BACKWARDS  INTEGRATION  LIMIT.  IF 

NO  BACK  INTEGRATION  IS  'EOUIRED  SET  T0 

TO  A  VALUE  GREATER  TRAN  OR  EQUAL  TO  Tl. 

TC.  —  COMMENCEMENT  TIME  FOR  THE  FILTER. 

U.T.  ME  AS > RED  IN  SECONOS. 

TE,  ENDING  TIME  FOR  FILTERING. 

U.T.  MEASURED  IN  SECONOS. 

TF,  FORHARO  INTEGRATION  TIME  LIMIT.  IF 

NO  INTEGRATION  IS  PF-jijjogo  JET  TF 
TC- A.  VALUE-  LESS  THAN  OR  EQUAL  TQ  T  E. 

TG,  DATA  REPROOUCING  TIME  LIMIT.  IF 
NO  REPROOUCING  IS  REUTREO  SET  TG 

TO  A  VALUE  LESS  THAN  OR  EQUAL  TO  TF. 

TX  TIME  SPAN  FROM  TC  ON  FOR  IMPROVING  LAUNCH 

CONDITIONS  WHEN  KEY39T*l  <  3EFAULT*  5  SEC.  } 


Figure  4b.  Card  Input  Format  for  DRIVEB 
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FILTERING,  INTEGRATION,  REPRODUCING  OPTION  KFV.  T* 

IS  KEY OPT  UNE (11  TO  PRO VI OS  NORMAL  TNIEG . -FT IT.  “UN, 

NINUS  ONE(-l)  FO»  ONE  PASS  INT  EG.-  FILT .  RUN, 
ZEROTH  TO  PROVIDE  [ N  T  E5R  ATIQN  WIT  POUT  DATA, 
TWO<2)  TO  PROVIQE  OATA  REPRODUCING  RUN, 

THREE (31  TO  PROVIDE  SINGLE  FILTERING  RUN, 

FOUR<4)  TO  PROVIDE  ODJPlE  FILTERING  RUN, 

USING  INlflAL  CONDITIONS  AT  TINE  Tl. 

INITIAL  CONDITIONS.  3C20.8 

17-  PIU-3J-  LN1TIAU  CONDITIONS  IV  GEOCENTRIC 

COORDINATES  AT  TINE  Tl. 

ALL  VALUES  IV  KILOMETERS. 

OR,  IF  »l  (1)  *P1<  ?)  =  »1  m  =C.  , 

aaiAiN  initial  conditions  at  launch 

FROM  LAUNCHER  LOCATION  OATA. 

18-  Pl<4-6)-.  -CORRESPONDING  VELOCIT.r  CONDITIONS* 

ALL  VALUES  IN  KILOMETERS /SECONO. 

OR-  ONLY  IF  Pl<6>*0.» 

AZCL,  AZIMUTH  of  ROCKET  CENTE’-LINE  AT  L  A'JNFH  OEG 

ELCL  ELEVATION  QF  ROCKET  CENTER-LINE  AT  LAUNCH  OEG 

NOISE  COVARTANCE  DATA.  fli.t 

19  DFACT  MULTIPLICATION  FACTOi  FOR  ASSUMED 

DRIVING  NOISE  COVARTANCE. 

PRINTEO  OUTPUT  CONTROL  CARD.  eT5 

20  IUNIT,  <-1  FOR  FEET  AND  f££T/SEC. 

♦  2-  FOR  KILOMETERS  AND  KM. /SEC. 

ISKIP,  INCREMENT  AT  WHICH  DATA  POINTS 

-  ARF  TO  at  PRINTED  AvO  PLOTTED*. 

rsAve,  selects  the  velocities  amo  stanoaro 

OEVIATI ONi  TO  9E  PRINTED. 

*1  FOR  GEOCENTRIC  VELOCITIES  ANO  hEVtattonS. 

♦  2  FOR.  VELOCITIES  REr  ER^NCEO  TO  THE  LAUNCHER. 

*3  FOR  VE.OCITTES  RffcR£ NCED  TO  THE  RAQAR. 

TCOM,  - COMMON- OHS P  SKIP  PARAMETER. 

POSITIVE  TO  GET  H9°SVI4T£0  DUMP, 

Nt-GATI V E  ro  GET  EXTENDED  DU-P. 

IOPT  TAPER  WRITE  OPTION  SE.ECTION 

a.  OR-3LANK  FDR  EVERT  ISKIP  POINT 
1  FOR  EVERT  POINT 

LAUNCHER  LOCATION. —  .  -  2C10.4 

21  AALTM  LAUNCHER  ALTITUDE  APDVE  THE 

SPHEROID  IN  METERS. 

22  ALOQEG,  LONGITUDE  IN  DEGREES,  MINUTES 

ALOKIN*  ANC  SECONDS.  SAME  CONVENTIONS 

AUDSEC  AS  RADAR  .ONGITUOE, 

23  ALAQEGv  .  ..LATITUDE  IN..  (DEGREES,  MINUTES 

ALA  MIN ,  ANO  SECONDS.  SAME  CONVENTIONS 

ALASEC  AS  RADAR  .ATITIUOE. 

PL" TTED  OUTPUT  CONTROL  CAROS.  T2,SX,7A10 

24  NVAR,  NVAR  SELECTS  THE  VARIABLE  TO  =E 

TTITLE(4I  PLOTTED  VS  TIME  frq«  LAUNCH  FPOM 

..THE. CON  TEXTS  OF  TAPE4.  THE  VALUES 
ALLOWEO  FOR  NVAR  CORRES°ONO  TO  THE  HOF" 

.  LOCATIONS  GIVEN  9 SLOW  FOP  TAPE 4. 

YTITLE  IS  THE  OPDINATE  * ITLE  ANO  t' 

LIMLTED  TD  40  CHARACTERS  . 

25  NG,  NG  SPECIFIES  THE  HUNTER  OF  CHARACTERS  CF 

3TI TLE  <  7)  G  TITLE  TO  8E  PLOTTED.  GJ  ITLE  hay 

CONTAIN  ,'JP  TO  70  CHARACTERS. 

28  ON  SELECTS  OTHER  VAPIA9UES  TO  9E  PLOT  TFP 

. —  EACH- WITH  NVAR^TTITL;  A»  FOR  CARD  ?4. 

NOTE.  ONLY  ONE  CARD  OF  TYP;  25  IS  ’E0JI°iD.  THE  P. mT  SEOuEmCE 
IS  AUTOHmTICALLY  TERMINATED  WHEN  THE  FOF  IS  e SCOUNTERED . 


FIGURE  4c.  CARO  INPUT  FOR  DRIVEB 
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Table  1.  Rockets  Modeled  for  DRIVEB. 


r 


Rocket 

Key 

AEROBEE  150 

I 

03 

AEROBEE  170 

04 

NIRO 

07 

NIKE  TOMAHAWK 

08 

UTE  TOMAHAWK 

09 

PAIUTE  TOMAHAWK 

10 

NIKE  HYDAC 

11 

NIKE  ORION 

12 

BLACK  BRANT  IV-B 

16 

BLACK  BRANT  V  -A 

17 

j  BLACK  BRANT  V  -B,  C 

18 

JAVELIN 

19 

NIKE  JAVELIN 

20 

ARIES 

24 

SERGEANT-EXCEDE 

25 

ASTROBEE  D 

30 

AEROBEE  350 

35 

CASTOR 

40 

CASTOR  LANCE 

41 

TALOS  CASTOR 

51 

ARBITRARY  UNPOWERED 

NONE 
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Figure  5.  Description  of  Filter  Time  Limits,  KEVOPT,  Geocentric  Position  and  Velocity 

C  -l  A  ~  i  _  onrwrn  J 


conditions  are  required.  The  numerals  after  the  process  INTEG,  FILTER, 
REPROD,  or  VFILT  indicate  the  sequence  in  which  the  trajectory  is  computed. 
Following  the  end  of  the  main  segment  at  TE,  integration  and/or  reproducing 
may  be  selected  by  TF  and  TG.  In  all  reproducing  cases,  and  for  the  single 
filtering  case  from  TC  to  TE,  the  velocity  estimate  only  is  smoothed  in  the 
backward  pass.  The  final  segment  goes  backward  from  TC  to  the  back 
computation  limit  TB.  KEY0PT=3  provides  a  run  identical  to  the  original  DRIVE3 
run,  and  KEY0PT=4  matches  the  original  DRIVE5.  KEY0PT=1  is  normally  used  now. 
When  KEY0PT=1  or  -1  apparent  radar  data  is  generated  from  launch  to  TC,  where 
good  actual  radar  data  becomes  available.  The  subsequent  backward  filtering 
to  TLNCH  results  in  a  consistent  trajectory  from  the  launch  pad  and  through 
the  radar  data.  When  KEY0PT=1,  launch  conditions  THFACT,  DFACT,  AZCL  and  ELCL 
are  optimized  to  obtain  the  smoothest  connection  at  TC,  TC  +  TX/5  and  TC  +  TX 
seconds.  If  this  search  procedure  is  unsatisfactory  for  some  reason  (such  as 
negative  thrust  delay),  it  may  be  necessary  to  select  a  different  connection 
time  TC,  or  to  fix  launch  conditions  and  inhibit  optimization  by  setting 
KEY0PT=-1 .  Notice  that  for  KEY0PT=1  or  -1  and  Pl(6)=0,  PI (4)  and  PI (5)  are 
taken  to  be  the  launch  azimuth  and  elevation  respectively  in  degrees,  further, 
PI ( 1 -3 )  are  ignored  and  the  launch  coordinates  are  derived  by  DRIVEB. 

Model  Uncertainty  -  QFACT 

QFACT  offers  a  general  method  for  controlling  the  significance  of  the  radar 
data  relative  to  the  model.  The  usual  range  for  QFACT  is  0.01  to  0.1,  with 
the  higher  values  indicating  greater  uncertainty  in  the  model  and  therefore  a 
greater  weight  on  the  data. 


Diagnostics 

A  diagnostic  print-out  is  provided  every  ICOMth  computational  step,  typically 
once  every  10  seconds  by  setting  ICOM=100  for  a  0.1  sec  data  rate,  and  at  a 
few  points  about  every  transition  stage  such  as  TC,  TE,  etc.  The  abbreviated 
dump  includes  geocentric,  radar  and  launcher  referenced  position,  velocity, 
and  acceleration  data  in  cartesian  or  polar  coordinates,  as  well  as  ballistic 
and  geographic  data.  The  extended  dump  with  negative  ICOM  also  includes 
various  transformation  and  covariance  matrices. 
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Output  TAPE4 


Complete  data  for  recreating  printed  and  plotted  reports  and  for  further 
statistical  computing,  if  required,  are  written  in  104-word  header  and 
data  binary  records  on  TAPE4.  The  format  and  contents  are  described  in 
Figure  6.  Launcher  referenced  coordinates,  ballistic  and  local  velocity 
data,  and  observed-computed  residuals  provide  useful  additional  information 
that  was  not  available  originally.!  TAPE4's  resulting  from  various  radar 
DRIVEB  solutions  of  a  rocket  trajectory  may  be  combined  two  at  a  time  in 
DRIVEC,  each  pass  producing  a  multi-radar  TAPE4  trajectory  estimate.  The 
printed  summary  report  is  designed  for  8-l/2"xll“  binding,  and  includes  a 
variety  of  information  contained  on  TAPE4  under  partial  control  of  ISAVE. 

2.3  DRIVEC 


Two  TAPE4  solutions  of  the  same  trajectory  are  input  to  DRIVEC  as  TAPE10 
and  TAPE11.  The  reference  trajectory  for  header  and  residual  purposes  is 
TAPE10.  Information  on  the  header  is  extensive,  so  that  only  minimal  print 
and  plot  control  data  needs  to  be  entered  on  cards  (Figure  7)  to  initiate 
the  multi-radar  solution.  Note  that  ISKIP  here  refers  to  the  increment 
TSKIP  at  which  TAPE10  was  generated,  and  not  to  the  original  computational 
step  DT.  There  is  presently  no  provision  for  unequal  TSKIP  increments 
on  TAPE10  and  TAPE11.  Another  item  to  note  is  that  DRIVEC  computes  the 
position-velocity  state  vector  by  statistical  weighting,  and  not  according 
to  any  trajectory  dynamics.  Therefore,  discontinuities  are  possible  where 
TAPE10  or  TAPE11  filtering  modes  change,  and  these  can  be  improved  by  revised 
DRIVEB  runs  with  modified  statistical  parameters  (see  Appendices  C,  F  and  G). 
The  final  TAPE  4  and  printed  and  plotted  output  are  identical  to  DRIVEB. 

2.4  REGEN 


The  printed  and  plowed  output  of  DRIVEB  or  DRIVEC  can  be  fully  reproduced 
using  TAPE4  in  prog; am  REGEN.  Header  information  is  comprehensive  so  that 
the  print  and  plot  report  generation  is  accomplished  with  minimal  card  data 
(Figure  8)  that  is  similar  to  the  DRIVEC  input.  ISKIP  here  refers  to  the 
original  or  apparent  computational  step  DT. 


OUTPUT  REPORT  OAT  A  FILE  NRITTEN  ON  TA®£4.  THE  FIRST  RkCppO 
CONTAINS  FLIGHT  INFORMATION  FOR  use  TV  PROGRAM  RE3FN,  mmTCH 
REPRODUCES  PRINTED  ANO  ®LOITeO  OUT’ur.  Tm£  REMAINING  RECORDS 
contain  data  to  bf  printed  and  plotted  ®v  either  ppjvt®  routines 

OR  RV  RSGEN.  IF  IOPT»l,ALL  G  AT  A  POINTS  ARE  NRITTEN  rN  T»PE4, 

IF  I o» T a 0 •  EVERY  ISKIPIH  point  is  NRITTEN. 

FIRST  RE'ORD 

MORO(S)  SYMBOL  OESCRIPTION 

I- 6  TITlEd-eidC  CHARACTER  ALPHANUMERIC  IOENT  IF  ]CA' ION 

9-10  OAT  (1-21  RUN  OATE,  9LANK 

II- 14  T  ITLE2 ( 1-4FURTHE R  ALPHANUMERIC  ID ENTIF1CATI  CN 

15-16  r ITLE2 ( 5-6PAYL0A 0  (KG),  RAREA  IKM»*2) 

IT  0  AT  ( 31  -CO 

18  DAT  ( L>  LAUNCH  OATE  (MO/DA/vr) 

19-20  DAT  (5-6)  ROCKET  ID 

21-21  OR0GI0(l-3PR0GRAM  ID  FOR  PLOTTED  OUTPUT 

24  HRAO  RAOAR  ALT ITUOE  (KM) 

25  RLOOEG  RAOAR  LQNGITUOE  IOE3) 

26  R  LON  IN  RAOAR  LONGITUQE  (MT'II 

27  RLDSEC  RAOAR  LONGITUDE  (SET) 

28  RLAOEG  RAOAR  LATITUDE  (DEG) 

29  RLANIjA-  RAOAR  LATITUDE  (MIN) 

30  RlASEC  RAOAR  LATITUDE  (SECT 

31  A ALT  LAUNCHER  AlTITUOE  <<N) 

32  alooeg  -mtncher  longitude  nr-i 

11  ALOMIN  LAUNCHER  LONGITUDE  (MIN) 

34  ALOSEC  LAUNCHER  LONGITUDE  (SEC) 

35  ALAOEG  .  LAUNCHER  LATITUDE  (DEG) 

36  ALAMIN  LAUNCHER  LATITUDE  (IIP) 

37  AlASEC  LAUNCHER  LATITUDE  (SEC) 

is  tench  launch  time  ut  (secs) 

39  THRAT  THRUST  RATIO-  to  nominal 

lc  delth  delay  in  thrust  (secsi 

41  ATCL  -- ROCKET- LAUNCH  AZIMUTH  (OEGS'T  .  - 

42  ELCL  ROCKET  LAUNCH  ELEVATION  (OEGS) 

43  31  DONPUTATLON  TIME  SIR®  SIZE  ISECS) 

44  TlI  INITIAL  PRINT  THE  ®RCM  LAUNCH 

45  TLF  FINAL.  PRINT  TIHE  FROH  LAUNCH 

44  TSKIP  TAPE  NRITTEN  EVERY  TSKIP  SECS. 

■♦7-104-  ZEROES-  ZERO-FILLED. 


S  US  SEQUENT  RECORDS. 

NPR  p  (3)  SYMBOL  DESCRIPTION 

1  GMT  UNIVERSA.  TINE  MOURa) 

2  VINE  TIHE  AFTER  LAUNCH  (SECONDS) 

3  SN.  iIGNAC  TO  NOIaE  RATIO 

4  RAG  RIGHT  ASCENSION  OF  SREENNICH  (RID  I 

5-7  ®2  ( 1-3)  FIUT.  GEOCENTRIC  PDStTION  /ECTOR  (K») 

O-l"  ®2 ( 4-6)  FILT.  GEOCENTRIC  VE.DCITY  VECTOP  (Km/«c> 

11-13  ®2  (7-9)  FILT-  GEOCENTRIC  ACCEL.  VECTOR  (K*/SEC»*2) 

14-16  ® VL  ( 1-3)  FILT.  lAJNCHER  REF.  POSITION  (KM) 

17-ID  PVL14-6)  FILT  LAUNCHER  REF.  VELOCITY  (KH/  EEr  > 

20-22  ° VL  ( 7-9 )  FILT.  lAJNCHER  REE.  ACCEL.  «H/«-«»2) 

23-25  3VL(1-3I  FILT-  LAJNCHER  REF-  RANGE, AZ, EL  («»,Cfc> 

26-21  DVL  ( 4-6)  FIlT.  LAJNCHFR  RRP.  R, A, e  RATES  (K“/*\  DEG/S) 

29-64  XSK 1-361  COMPLETE  ERROR  COVARIANCE  MATRIX 

65-67  0  VR  ( 1-  3)  FILT.  RA3AR  °ANG”.AZ,  EL  (KM,  PEG) 

68-70  DVR  (4-6)  -  F.ILT--  HAJAP  RANGE -AZ,  EL  RATES  (  KM  /S.OEG/S  ) 

7 1  - T3  0 VR  ( 1-3)  FILT.  RAOAR  REF,  POSITION  VECTOR  (KM) 

74-75  S P V ( 1- 3 )  GEODETIC  ALT.,  LONG(H),  LAT.  (KM,  CEGI 

77  VR  MAG.  OF  LOCAL  VELOCITY  VECTOR  (Km/SEP) 

78-79  AZS.ELR  A7, EL  OF  LOCAL  VELOCITY  VECTOR  (D  EG) 

80  CO  DRAG  COEFFICIENT 

8 1  DEN  ATMOSPHERIC  DENSITY  (KG/KM**3> 

82-63  ORAG.TH  DPAG,  THRUST  FJR'-ES  (KDKH/S2) 

84-85  DVR  ( 1-3 )  RAN  RAOAR  DATA  RANGE  ,  A  Z ,  EL  (KM, DEG) 

8  7-  ®  9  D  VG  (  1-3)  GEOCENTRIC  PAN  RAOA®  DATA  VECTOR  (<“> 

90  GRANGE  GROUNO  RANGE  ALONG  SPHEROID  (KM) 

91  1CLMAG  LAUNCHER  REE.  ACCEL.  “AGNITUOE  (K7V«2) 

92  ®c.SALl  ALTITUOE  “ESIDUAL  (<M)  .  _ 

91-95  <RSG,YRSG,ZRSG 

GEOCENTRIC  »OSIIION  VICTOR  FESIDUAL'  (KM) 

9®  DRSG  GEOCENTRIC  VECTOR  MAGNITUDE  RESIDUAL  (KM) 

97-9  9  vRSL.YRSl.ZPSL 

launcher  position  vector  residuals  (k*m 
108  GRSL  LAUNCHER  VECTO®  MAGNITUDE  RESIDUAL  (K»> 

L 0 1 -ID T  KRSR.YPSR.IPaR 

RADA®  POSITION  V  ECT  D  ®  NSSIDUAlS  (  »“) 

104  GRSR  ■ RAOn®  VEDTOR  MAGNII  ))F  RESIOuAL  (KM) 

NOTEl  all  PESirUALS  APE  (PAN-EILI ERSD) 


Figure  6*  TAPE4  Header  and  Data  Records--Output  by  DRIVEB  or  DRIVEC 

Input  for  DRIVEC  or  REGEN 
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OROG-.AN  ORIVECI I  NF  UT  , OUTPUT,  T  A  f  E  A  , TAPE7,TAPEi;,TAPEii, 
t - -W-6 -WT c+H  (  T  AfEF^TH  W  ) - 


-iwe: - &€vst-0PS,?-i?7— 

- OTrr»tPTTCFr~- - 


version  -sep-.  ,1573 — 


I/O  "Efihit ion; . 


CATS  INPUT  FCOCUCEC  AS  OUTPUT  FROM  F,Ivgo 

PRINTED  AND  CLQTTjr  OUTPUT  TAT  A  CONTAINS 
3gtgCT?fr  R-eihT-?"  IF1  ;0PT*".-'  ALL  PfrlNTT  "I»~IO TH 


ITQRAPS 


TAFE7  SC'ATCH  FJLE  ECR  DATA  M  ANT  PULA  T  ICN .  NASS 

- ACCESS  ORGANIZATION.- - - - 

TAPE*  SYSTFM  PRINTER. 

Fg-AD-gR  ■. 


0  ATA—TNPDT'  CS^MPIPT!;>r 


—5 - rt - M - - - ** - 

CAPO  ITF»  DESCRIPTION  F  CPN  A T 


r-A RO  ♦ - PRINT  AN?  0eTPi'T--C-&NTC0L---=ATAT - — - 5-1= - 

IUNIT,  *1  FOR  F£ET  AND  fEET/SlCONN 

- -  rop  F-I't'tmgTFRS— ANf>— SH  ./=■£»■. - 

I^XIP,  INCREMENT  at  WHICH  data  POINTS 

- ARC"  TO  eg-  PRINTED  ANO-  PtOT^D-. - - -  - 

ISAVE,  SELECTS  THE  VELOCITIES  AND  ST  AN  C  ARO 

- DEVIATIONS  to  BE- -PRINTED. - — - - 

♦1  FOR  GEOCENTRIC  VELOCITIES  ANC  DEVIATIONS. 

- +9-+Q*-  VELOCITIES  R£  FE'96  NO&O'  TO  ■  TH=--L  AUwCHFC  . - 

♦3  FOR  VELOCITIES  REFERENCED  TD  THE  RADAR. 

- 1  EG*, - SONHON-  8-UHP-  GW  I  F— &AR  A-HET  ER  r-SER  yg~.-c  - 

FUNCTION  AS  TSKIP  PARA^ETFC. 

- fAPT - T  A  PER  LP1T5-0FT  I6N--5ELE&TI6-N - 

0  OP  ELANK  FCS  EVERY  ISKIP  POINT 

- 1  .  r-Qft  ■  CV=°Y-  POINT - 

OAPP--J - NAHEVIST  -OAT A  . - T-L  ItTL-F,  -*M€-  WIN-T-T-INE-L-W-Ts-HAS— »  f-F— 

AFTER  LAUNCH)  HAY  9E  fESET.  FOR  EXA-cLi, 

- p - t - 3 - 3 - a - 5 - 

111*. 

- SNAH'A  Tt-I«-tDvO  50 .  OS - 

IF  NOT  RESET,  rtf  ORIGINAL  LIMITS  HILL  3£  USE1". 


TITLE  INFORMATION  IS  INFUT  from  the  HEADER  RECORD  ON  tafgio 
-HHI-CH-NAS-  the  -T-APEV-^UTPUT-  PRCN-  PRI  ve°. - - - — 

■FEORT-9— WTPUT- CONTROL-  CAROS. — FCP»»T  -IE,  3X  ,  *«  t3 - 

I  NVAP,  SELECTS  the  variable  r0  36  ’LOTTED  FROM  thc 

- Y-T-I-TLE-*  1— -A  )  CONTENTS  0=  TAPgfc.,  YTITLE  afccA  c  S-4S— FH€ - - 

ORDINATE  TITLE.  (LIHIt  OF  U 0  CHARACTERS) 

- A — w-G, - saJS-IRfE-S-THE-HUH^ER  OF— SWWA^TE-55— I N-  THE — --  - 

Gr:TLE(l-7)  GENERAL  title.  GTITLE  APPEARS  ON  EACH  CL'TPUT 

- 6RAFH,  <ha»IHUW'  6F  ■? 5--CH AR A6Ti£-R G ) - 

CARP  5  ON  SELECTS  other  VARIABLES  TO  3 E  FLCTTEP 

- - - EACH  WITH->WAR, ■YTITLg  AS  CDR-CACG— Zr - 

NOTE.  Only  ONE  capo  OF  TYPE  ^  IS  PECUIPEO.  the  PLOT  SEOUENPS 
IS  AUTCHATICALLT  tcrhINATEC  WHEN  THE  EOF  IS  ENCOUNTERS. 


-gU-ictfr—c  gaopr  DATA  FILE  wc-ITTEN-CN-  TAP^A^-FHE- «fis-F-T— oroCR= - 

CCNT A INr  flight  INFORMATION  FOR  USE  3Y  proGRA-  REGEN,  WHICH 

-  Brae OQ u*g»'  PRINT EO  AND  ct^T^O  OUTPUT t  TuC  RE M » ! NlN«-  RF  OCR  09 - 

CONTAIN  DATA  ’c  ?E  FRINTEO  JNO  PLOTTED  EITwf?  DRIVEC  ROUTINES 

IF  IOPT=0,  EVERY  ISKIPTM  POINT  is  written. 


Figure  7.  Card  Input  Format  for  DRIVEC 


-23- 


c FOG? AM  PE CEMTNFUT,  OUTPUT  ,T  APE’S, 
*»*.*  *»**■» 

P  -  GEN  0  EVSLOPEO  RV  LCGICGN, I NC . 


TAPE 4,T APES- INPUT, 1 AF  E5=  OUTPUT) 
VERSION  SrP.,l973 


-■z^e^rf  S^f^-WK-eR-r-  G  E  NS-^A"TE‘S'-eitY*,t}Jf — TRA  JrC-TOfrV-  < ■•-°DC-TS — 
USING  T  A  c£4  .  INPUT  DATA  CFODUr=;n  3V  H  c  I  VE  3 


THE  FOLLOWING  CARD  INPUT  IS  EEOUIFEH 

CARD  1-  PRINT  CON*Ret~OATA -  —  FORMAT-41-  .  . 

IIIN3T,  +1  FOR  FEET  AND  FEET/SECONO 

IEKIP,  DATA  FCINTS  WILL  RE  PRINTED  EVERY  ISKIF*CT  SECONDS 

TCAV5 — SEL£CT-S  TME--VELC  CITIES  A  NO — RE-SI  C-UAtS  T  G— EE— P-RINTE0 - 

♦ 1  FOR  GEOCENTRIC 

■  ---  . — ►  2- FOR- t-MiHe+«R  - - -  '  -  -  - 

+3  FOR  RADAR 


CARO  2  NAMELIST  OATA,  TLI,TLF,  the  PRINT  TIME  LIMITS ( AS  TIME 

AFTER  lAUNOH)  MA-Y-BE--RESET .  FOR  EXAMPLE-*  - . 

0  12  3  4  5 

1  1  . . t~  i  -  1  1 

SNAMA  TLI=15.0  ,TLF=450.  0? 

- IF  NOT  ReSE-Tj-TT-t-O-R^-GlNAL  LIMlTS~-HiEfatr- >r  OS-D? - 

TITLE  INFCR NATION— I S  INPUT  FROM  THE  WEAOE  R— RECORD  ON  TAFE4-. - 


PLOTTED  OUTPUT  CONTROL-EA-PRS. - F CPMA-T — 12,  8X  ,  7  e  !•]  - 

3  NVAR,  SELECTS  THE  VARIABLE  TO  BE  PLOTTED  FROM  THE 

YT-T-T-tc-Hr»4->— CTH^Tg-NTS-  OP — TA  Pr  4  .  YTE-Tfc-g  -APPEALS-  AS — T HE - 

OROINATE  TITLE.  C L IMIT  OF  40  CHARACTERS) 

i,  N&, _ _ £pKGIF-J£5-THP  NUMBER-  OF  CHARACTERS  IN-  THE - 

GTITtE  (1-7)  GENERAL  TITLE.  CTITLE  APPEARS  ON  EACH  OUTPUT 

.. _  _ G^Am.  (HA  > I-MU H-  -  C  F-  7 0— CH-ARACTERS) - - - 

card  5  ON  SELECTS  OTHER  V ARI £3LES  TO  BE  PLOTTED 

NCTE.  RFFER  TO  THE  USER'S  MANUAL  CR  THE  ORIYER  LISTING 
^F^qihfR  0 ETAEL-S- G-N- -THE— SETU p  FOR-TPE SE- L- AS-T- Pfc Q-T- 
ALL  FILF  DEFINITICNS  ARE  GIVEN  IN  THE  ORIVE9  LISTING. 


FCR 

-CAPOS. 


Figure  8.  Card  Input  Format  for  REGEN 
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2.5 


Operational  System 

The  rocket  trajectory  operating  system  has  been  modified  to  reduce  the 
need  for  complete  binary  programs  of  DRIVEA,  DRIVES,  DRIVEC,  and  REGEN. 
The  SCOPE-EDITLIB  facility  is  used  to  provide  one  binary  library  for  all 
subroutines.  The  main  program  with  the  LIBRARY  control  card  then  causes 
unsatisfied  externals  to  be  loaded.  Program  maintenance  thus  requires 
only  one  file  of  source  routines  in  UPDATE  form. 
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3. 


PROGRAM  GUIDE 


in  this  section  the  organization  and  program  flow  of  the  four  separate 
programs  are  described.  The  subroutines  that  are  called  are  generally 
modular  and  self-explanatory.  Most  of  these  routines  are  described  in 
reference  1;  revisions  and  new  routines  are  described  below  and  in  the 
appendices. 

3.1  DRIVEA 


This  is  the  raw  radar  data  editing-preprocessing  program.  The  original 
system  used  a  separate  statistical  preprocessor  DRIVE2  to  clean  up  data 
values  for  the  next  stage  in  DRIVEB,  and  prepare  initial  plots  of  the 
trajectory.  The  estimator-filter  proved  to  be  an  excessive  user  of  CP  time 
and  memory.  Further,  the  predictive  filter  tended  to  ring,  giving  poor 
position  and  velocity  estimates.  All  editor  and  preprocessor  functions 
have  been  combined  in  DRIVEA.  A  five-point  running  average  of  accepted 
data  serves  as  the  velocity  estimator.  A  complete  report  on  the  editing 
and  estimating  is  printed.  This  program  produces  results  that  are 
consistently  superior  to  the  original  DRIVE1  and  DRIVE2  while  requiring 
less  memory  and  about  40  percent  of  the  computer  time.  Additional  features 
have  also  been  added  to  facilitate  processing  tapes  with  unreadable  data 
segments  by  the  use  of  spanning  polynomials.  Special  buffering  and 
mass-storage  to  tape  file  conversions  which  were  carried  out  with  CPRS  are 
unneccessary  and  have  been  eliminated. 

Pre-edit  Scan 

Data  on  cards  in  the  format  described  in  Section  2.1  are  read  and  printed. 
Pages  111-3,4  of  Vol.  I  of  Reference  1  are  reproduced  in  Figure  1,  and 
describe  the  raw  input  TAPE1{=RAWDAT)  header  and  data  records.  NSKIP  6-sample 
records  are  skipped  and  the  (5-N0K)th  and  (6-N0K)th  samples,  i.e.,  for  default 
N0K=0,  the  5th  and  6th  samples,  are  taken  as  the  base  starting  points  for 
initial  position  and  rate.  A  complete  scan  of  the  input  tape  is  first  made  to 
determine  all  quadrants  that  the  azimuth  covers.  Logic  is  implemented  to 
provide  numerical  continuity  if  the  azimuth  goes  through  zero,  i.e.,  360° 
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jump.  In  the  unlikely  event  that  the  azimuth  crosses  both  0°  and  180°,  the 
present  logic  is  inadequate  and  this  is  flagged  as  requiring  a  special  fix. 

Editing  and  Spanning 

Following  this  pre-edit  scan  the  tape  is  repositioned  and  the  main  editor 
pass  is  carried  out  to  create  a  clean  file,  free  of  gross  or  unreasonable 
sample  values.  The  operation  of  the  editor  is  described  in  detail  in  Appendix 
A.  The  cleaned  file  and  associated  flags  for  corrections  that  may  have  been 
made  are  saved  on  temporary  sequential  file  TAPE11.  At  this  point  a  special 
pass  is  to  be  made  to  determine  smooth  spanning  cubic  fits  over  individually 
selectable  ranges.  These  curves  are  introduced  in  the  appropriate  segments 
during  the  final  preprocessing  phase.  (See  Appendix  A. 3  for  details) 

Output  and  Geocentric  Estimates 

Creation  of  the  edited  output  file  TAPE10  is  carried  out  over  the  range 
TLO-THI.  Header  and  data  format  are  described  in  Figure  3.  Raw  geocentric 
velocity  and  position  are  estimated  at  the  same  time,  are  written  on  temporary 
files  TAPE8  and  TAPE9  respectively,  and  are  printed  out  along  with  the  edited 
range,  azimuth,  elevation  values.  Subroutine  FINDER  calculates  the  geocentric 
position  corresponding  to  the  radar  data,  and  also  obtains  local  geodetic 
coordinates.  Altitude  is  used  to  determine  refraction  corrections,  if 
required.  A  simple  velocity  difference  clipper  and  5-point  velocity  smoother 
is  used  to  estimate  the  geocentric  velocity  vector.  This  results  in  an  offset 
of  2  samples  which  is  not  serious  but  is  corrected  for  by  a  hold  buffer.  The 
geocentric  state  vector  is  required  by  DRIVEB  as  the  starting  estimate  for 
filtering  from  any  time  after  launch.  The  print-out  for  each  sample  includes 
time,  edited  range,  azimuth,  elevation  with  correction  flags  if  applicable, 
and  from  TLO  to  THI,  altitude  and  the  estimated  geocentric  position  and 
velocity  vectors.  Finally,  subroutine  ARTIST  provides  a  plot  of  these 
variables  versus  time,  using  Tapes  8,  9  and  10. 

3.2  DRIVEB 


This  is  the  main  filtering  and  trajectory  generation  program.  Data  on 
cards  in  the  format  described  in  Section  2.2  are  read  except  for  the  plot 
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control  cards  which  are  read  after  completion  of  processing  and  filtering. 
The  header  of  the  TAPE10  edited  radar  data  provides  creation  data  for  this 
filej  this  overrides  card  input  data  except  when  purely  integrating  (KEY0PT=0). 
The  maximum  capabilities  of  DRIVEB  are  exercised  for  the  launch  optimizing, 
filtering,  and  terminal  reproducing  run. 

Time  -  Tape  Registration 

A  one-to-one  correspondence  is  first  established  between  sample  times  on 
TAPE10  and  internal  TAPE7  random  access  record  numbers  (MT).  Requested 
time  ranges  for  filtering,  integrating,  etc.  are  re-established  in  terms 
of  the  actual  sample  times.  The  input  tape  is  positioned  after  reading 
the  sample  for  the  initial  or  the  first  computation  time,  whichever  is 
greater,  and  is  read  again  as  the  processing  is  incremented  to  successive 
sample  times. 


Derivation  of  TAPE7  Internal  Data  Base 

Parameters  to  complete  a  TAPE7  working  record  for  any  sample  time  are 
described  in  Figure  9.  The  derivation  of  these  parameters  which  comprise 
the  tota’  information  base  is  now  discussed.  The  signal -to- noise  ratio  SN 
is  obtained  from  RADAR.  If  pure  integrating  is  being  employed,  SN  is  set 
to  -100.  If  pure  reproducing  is  used,  SN  is  set  to  +100.  These  special 
values  of  SN  control  multi -radar  combination  modes  in  DRIVEC. 

Statistical  Filtering  in  the  Geocentric  Domain 

The  geocentric  state  posit ion- velocity  vector  P2 ( 1-6 )  is  the  Kalman  weighted 
resultant  of  the  integrated  dynamics  of  RK  and  the  observed  raw  geocentric 
state  vector  DVG(l-6)  on  the  forward  filtering  pass.  The  observed  position 
DVG ( 1 -3 )  is  calculated  in  EULER  from  the  raw  radar  range,  azimuth,  elevation, 
DVR ( 1-3 ) .  Co-azimuth  was  used  earlier  (Ref.  1),  but  has  been  eliminated.  The 
observed  raw  velocity  DVG(4-6)  is  estimated  in  EULER  as  the  difference  between 
position  samples,  limited  to  a  change  of  0.01  km/sec.  from  the  existing 
filtered  estimate.  The  raw  velocity  DVG(4-6)  tends  to  be  tolerance  limited, 
but  the  random  access  noise  is  smoothed  very  satisfactorily  by  the  subsequent 
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Figure  9.  Description  of  Internal  Working 
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filtering.  On  the  backward  filtering  pass,  maximum  likelihood  or  interval 
estimation  weighting  of  the  backward  integrated  dynamics  and  the  original 
forward  estimate  completes  the  computation  of  P2  in  DRIVEB.  The  error 
covariance  matrices  XSI  and  XSIP  are  weighted,  as  described  in  Appendices  B 
and  F,  and  calculated  in  KAL  during  the  forward  pass.  XSI  is  recomputed  in 
DRIVEB  during  the  backward  pass.  Geocentric  acceleration  P2 (7-9 )  is  estimated 
at  any  time  by  a  call  to  FOR  and  WPROD  using  a  nominal  time  increment  of  1 
second. 

Radar,  Launcher  and  Other  Parameters 

Filtered  radar  range,  azimuth,  elevation  and  rates  OVR ( 1 -6 )  are  derived 
from  P2 (1-6 )  by  TRACOR.  The  radar  referenced  position  vector  PVR(I-3)  is 
not  required,  but  is  calculated  later  by  GENER8  for  report  and  TAPE4  purposes. 
Filtered  launcher  range,  azimuth,  elevation  and  rates  OVL ( 1-6 ) ,  as  well  as 
launcher  referenced  position,  velocity  and  acceleration  vectors  PVL (1-9)  are 
obtained  from  P2(l-9)  by  a  call  to  ROTATE,  a  routine  which  has  been  extended 
to  provide  these  general  transformations.  Local  velocity  and  ballistic 
parameters  VR,  AZR,  ELR,  CD,  DEN,  DRAG,  TH  were  obtained  when  FOR  was  called. 
The  TAPE 7  working  record  is  designed  to  match  the  TAPE4  parameter  list  through 
word  89,  i.e.,  DVG(3).  The  completion  of  parameter  calculations  for  TAPE4  is 
done  in  GENER8. 

Processing  in  Different  Time  Zones 

TAPE 7  parameters  for  record  number  MT1 ,  where  the  filtering  or  integration 
is  started,  are  calculated  based  on  initial  conditions.  Estimates  for 
successive  times  are  based  on  a  combination  of  the  integration  estimate 
using  RK  (see  Appendix  D)  and  the  raw  radar  estimate,  as  discussed  above 
and  in  Appendix  B.  Filtering  is  performed  from  MT1  to  MT2,  forward 
integration  from  MT2  to  MT3,  and  reproducing  of  the  raw  radar  estimate  from 
MT3  to  MT4,  or  until  the  altitude  goes  negative.  Then,  interval  or  maximum 
likelihood  estimation  filtering  is  performed  backwards  to  MT1,  with  a  final 
backward  integration  if  requested.  Figure  5  shows  the  processing  sequence 
and  operating  modes  available  with  the  various  KEYOPT  options.  In  all 
reproducing  and  single  filtering  segments,  the  backward  filtering  VFILT 


smooths  the  velocity  estimate,  but  the  position  values  front  the  forward 
pass  are  retained.  Backward  filtering  of  integrated  results  from  MT2  to 
MT3  is  performed  if  subsequent  reproducing  is  required.  Except  for  VFILT, 
all  backward  filtering  stages  are  the  maximum  likelihood  kind. 

Launch  Optimizing  and  Filtering  Mode 

The  launch  optimizing  mode  (KEYOPT^l)  achieves  the  most  satisfactory 
trajectory  estimates  in  the  usual  case  when  good  radar  data  is  available 
only  some  time  after  launch.  An  iterative  run  is  first  made  to  optimize  on 
conditions,  viz,  thrust  ratio  to  nominal,  delay  in  onset  of  thrust  from 
nominal  launch  time,  launch  azimuth  and  elevation.  The  procedure  is 
implemented  in  routine  LEASQ  which  searches  for  the  best  position  and  velocity 
match  to  radar  data  near  time  TC  (see  Appendix  E).  In  most  applications  radar 
data  is  erratic  from  launch  until  after  burn-out,  since  rapid  slewing  and 
large  accelerations  are  encountered.  Therefore  TC  is  frequently  set 
satisfactorily  at  10  to  20  secs  after  thrust  burn-out.  Radar  data  before 
time  TC  is  ignored,  using  only  the  forward  integration  results.  On  the 
backward  pass  the  filtering  and  integration  results  are  uniformly  refiltered. 

Output  TAPE4 

DRIVEB  results  are  periodically  printed  out  in  abbreviated  or  extended 
form,  and  include  a  few  extra  frames  around  each  mode  transition  time. 
Following  completion  of  the  above  processing  under  any  of  the  KEYOPT  options, 
TAPE7  results  are  converted  to  report  and  storage  form  on  TAPE4  by  GENER8. 
WRITER  provides  the  printed  user  report  and  plots.  The  TAPE4  header  and  data 
contain  sufficient  information  to  create  such  reports  and  plots  at  some  time 
in  the  future,  using  program  REGEN  and  subroutine  WRITER,  with  a  minimum  of 
redundant  card  data  input. 

3.3  DRIVEC 


The  multi-radar  trajectory  program  DRIVEC  was  developed  as  a  revision  of 
program  DRIVE4.1  The  card  input  requirements  are  greatly  reduced,  using 
instead  the  header  information  on  TAPE4  generated  in  program  DRIVEB  (Figure 
6). 
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Scope  of  Program 


The  primary  purpose  for  the  multi-radar  program  is  to  combine  the  results 
of  two  separate  radar  observations  of  a  given  rocket  flight  to  give  a 
composite  trajectory.  The  mathematical  and  logical  procedures  that 
accomplish  the  combining  operation  on  to  TAPE7  internal  working  records 
are  discussed  in  Appendix  C.  The  methods  used  to  accomplish  this  data 
combination  are  linear;  therefore,  an  arbitrary  number  of  different  sets  of 
radar  data  pertaining  to  the  same  flight  may  be  combined  by  repeated  execution 
of  the  program.  The  graphical,  printed  report  and  file  outputs  of  this 
program  are  identical  in  all  respects  to  the  outputs  of  Program  DRIVES. 

Program  Functions 

Specific  functions  of  Program  DRIVEC  for  processing  two  TAPE4  format  input 
files,  TAPE 10  and  TAPE 11,  include  the  following: 

—Determine  the  type  of  data  region  being  read  for  each  file 

--Use  the  region  types  to  select  proper  mode  for  combination 
equations 

--Evaluate  combination  equations,  depending  on  the  particular 
mode,  to  produce  the  composite  observation  vector  and  co- 
variance  matrix 

--Complete  the  calculations  for  the  remainder  of  each  output 
file  record  using  composite  values 

—Generate  a  printed  summary  of  the  combination  process,  i.e., 
an  event  log  for  the  resulting  output  file 

—Produce  an  output  TAPE4  file  containing  the  composite 
trajectory 

--Generate  printed  and  graphical  reports  in  the  same  format  as 
DRIVES. 

For  the  purpose  of  calculating  ballistic  and  thrust  parameters  as  well  as 
residuals,  TAPE10  is  taken  as  the  reference  trajectory.  Thrust  ratio, 
delay,  etc.  are  obtained  from  the  header  and  used  in  a  call  to  FOR.  Raw 
components  of  the  radar  data  are  also  obtained  from  TAPE10  for  use  when 
TAPE4  is  completed  in  GENER8.  Note  that  I SKIP  and  IOPT  refer  to  the  increment 
TSKIP  at  which  TAPE10  was  written  rather  than  the  original  integration  time 
step  DT. 
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3.4 


REGEN 


REGEN  is  a  free-standing  report  generator  program  which  uses  a  previously 
created  TAPE4  and  provides  identical  reports  and  plots  to  those  obtained 
after  a  DRIVEB  or  DRIVEC  run.  The  card  input  requirements  are  greatly  reduced 
from  the  original  program,  using  instead  the  header  information  on  TAPE4, 
after  which  WRITER  is  directly  called.  As  in  DRIVEB,  1SKIP  refers  to  the 
integration  or  basic  time  step  DT.  The  procedure  recreates  conditions  after 
a  DRIVEB  or  DRIVEC  run  has  called  GENER8;  otherwise  REGEN  is  simply  a  subset 
of  those  two  programs. 


Appendix  A 


The  Editor 


The  DRIVEA  editor  processes  raw  radar  data,  which  may  contain  a  variety  of 
errors,  and  generates  a  clean  file  with  reasonable  range,  azimuth  and  eleva¬ 
tion  readings  for  each  sample  time  over  the  span  of  the  original  data.  Errors 
which  are  handled  Include: 

a)  Missing  time  points 

b)  Data  with  grossly  improper  or  extraneous  time  points 

c)  Data  with  marginally  improper  time  points 

d)  Grossly  improper  data 

e)  Marginally  improper  data  (based  on  wide  tolerances 
set  by  user) 

In  addition,  the  user  may  selectively  cause  spans  of  range,  azimuth,  or 
elevation  data  to  be  replaced  by  a  smooth  interpolation. 

A. 1  Operation  of  the  Editor 


Starting  with  a  pair  of  successive  observations  that  the  user  specifies  as 
acceptable  for  time,  range,  azimuth  and  elevation  values,  the  raw  radar 
data  is  checked  for  proper  time  values  and  reasonable  rates  of  change.  If 
an  observation  is  found  to  be  improper,  an  entry  is  started  in  an  edit 
buffer  with  the  intent  of  replacing  this  data  with  interpolated  values  when 
acceptable  data  are  subsequently  found.  The  edited  data  for  as  many  time 
points  as  were  adjusted  are  then  written  to  the  output  file,  the  edit  buffer 
is  cleared,  and  processing  continues  on  new  radar  data  with  the  latest  values 
and  rates  of  change  providing  the  acceptance  criteria. 

A  typical  sampling  rate  is  0.1  sec.  Hence,  if  the  next  sample  time  is 
within  0.05  sec.  it  is  considered  only  marginally  improper,  and  Is  accepted 
with  the  time  corrected  to  the  0.1  sec.  increment.  A  time  value  lower  than 
this  0.05  sec.  tolerance  or  unacceptably  high  (say,  6.75  secs,  or  more)  Is 
considered  extraneous,  and  the  complete  sample  is  discarded  as  suspect.  If 


-34- 


the  time  value  for  the  next  sample  Is  moderately  high,  observations  are 
probably  missing,  and  an  interpolation  is  carried  out  on  all  the  values 
provided  the  change  rates  are  proper  and  provided  the  next  raw  data  sample 
does  not  reveal  the  time  jump  to  have  been  due  to  an  improper  time  reading 
rather  than  the  results  of  missing  data. 

As  long  as  the  time  readings  are  acceptable  the  data  will  not  be  discarded; 
instead  the  range,  azimuth  and  elevation  are  individually  checked  against 
selected  tolerances  and  adjusted  if  necessary.  If  the  observation  differs 
from  the  predicted  value  by  more  than  the  tolerance,  the  data  will  be  replaced 
by  interpolation  to  subsequent  acceptable  data.  If  the  observation  is 
marginally  improper  as  determined  by  a  tolerance  factor  QT0L=0.25,  the  data  is 
accepted  after  making  a  compromise  between  the  value  and  the  prediction.  If 
seven  or  more  observations  have  been  rejected,  QTOL  is  increased  to  0.5  to 
avoid  missing  a  new  trend.  After  modification,  the  value  is  treated  as  an 
accepted  observation.  The  modification  formula  is  simply 

Modified  Obs  =  (1  -  QTOL)  *  Prediction  +  QTOL  *  Actual  Obs. 

This  editing  process  for  each  variable  is  accompanied  by  count  indicators 
which  serve  two  purposes:  (1)  if  non-zero,  it  gives  the  count  from  the 
previous  accepted  sample  for  interpolation  when  a  good  sample  is  again 
obtained;  and  (2)  the  print-out  of  the  final  output  file  uses  count 
indicator  status  to  flag  whether  a  value  was  added,  interpolated,  modified, 
or  unaltered,  (see  Figure  A) 

A. 2  Implementation  -  Example 

The  implementation  of  the  editor  is  described  below  with  the  help  of  an 
example  for  radar  data  containing  a  variety  of  errors. 

The  edit  buffer  is  maintained  for  radar  data  and  for  error  status  indication. 
A  zero  indicates  that  ro  change  was  made  in  the  data,  while  a  101  denotes  a 
minor  modification.  If  the  data  are  all  accepted  (row  only  has  0  or  101)  at 
any  sample  time,  it  shows  that  all  Interpolations  for  any  unacceptable  or 
missing  data  can  be  completed  for  previous  observations.  Then  the  buffers  may 
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Raw  Radar  Data  Error  Status  Counters  Edited  Radar  Data 
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Figure  A.  Example  of  Editor  Processing 
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be  written  to  the  output  file  and  cleared  for  continued  editing.  In  the 
example  in  Figure  A,  the  longest  retention  required  in  the  buffer  is  seen  to 
be  from  7.4  to  7.9  secs,  or  6  samples.  The  present  program  allows  up  to  100 
samples  or  10  secs,  of  editing  which  is  limited  only  by  the  size  of  the 
buffers. 

The  observations  are  seen  to  be  clean  at  7.1  and  7.2  secs.  Change  rates 
are  re-estimated  for  every  observation  with  a  10  sample  time  constant. 
The  predicted  values  are  compared  with  the  raw  values  using  tolerance  limits, 
usually  5  km  and  5°,  for  complete  rejection.  A  lower  tolerance  limit  (25%)  is 
used  to  fully  accept  the  data,  or  if  the  error  exceeds  this  amount  to  make 
compromise  modifications  as  discussed  earlier.  Illustrations  of  typical 
editing  actions  A  thru  H  follow.  (A)  In  the  example  this  type  of  modification 
is  seen  to  have  been  made  in  the  range,  azimuth  and  elevation  where  the 
status  counters  are  101.  (B)  The  sample  at  8.4?  secs,  is  judged  to  actually 

be  at  8.4  secs,  for  consistent  time  incrementing.  (C)  The  sample  at  7.92 
secs,  however  is  less  than  0.05  secs,  from  the  earlier  sample  and  is  discarded 
as  extraneous.  (D)  The  sample  after  7.4  secs,  has  an  excessive  time  jump  and 
is  completely  discarded  even  though  the  associated  observations  are 
reasonable.  These  values  will  be  interpolated.  (E)  The  jump  in  time  after 
7.2  secs,  and  7.6  secs,  also  call  for  interpolation.  This  is  done  in  both 
cases,  but  in  the  former  case  the  next  sample  at  7.4  secs,  shows  that  the 
interpolation  between  7.2  and  8.3  secs,  was  a  false  alarm.  This  is  replaced 
by  a  simple  interpolation  for  the  7.3  sec.  sample  only.  (F)  Grossly  improper 
data  such  as  the  elevation  values  at  8.3  and  8.4  secs,  must  be  replaced  by 
interpolation.  The  status  counters  are  correspondingly  incremented  from  1 
thru  99  (the  buffer  size  limit)  until  a  good  sample  is  encountered  (0  or 
101).  (G)  &  (H)  Since  this  interpolation  is  identical  to  the  case  for  missing 
data,  connected  instances  of  bad  data  and  missing  data  cause  the  appropriate 
counters  to  continue  incrementing. 

The  modification  procedure  bears  some  further  discussion.  A  slow  adjustment 
band  is  necessary  between  the  limits  of  complete  rejection  and  complete 
acceptance  so  as  to  be  influenced  by  possible  new  trends  without  being 
permanently  and  incorrectly  thrown  off  by  a  few  offset  samples.  In  fact, 
if  7  samples  in  succession  are  completely  rejected,  the  modification  tolerance 
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factor  is  increased  from  0.25  to  0.5  so  as  to  react  faster  to  diverging 
trends.  The  example  shows  the  rapid  adjustment  of  the  range  at  8.2  and  8.3 
secs,  to  allow  for  the  8.51  km  readings.  But  the  subsequent  value  of  7.13  km 
at  8.4  secs,  is  still  fully  acceptable.  The  wide  variation  has  been  greatly 
reduced,  but  its  smoothed  elimination  is  possible  only  with  the  subsequent 
filtering  program  DRIVEB.  Some  improvements  may  be  made  however  by  resetting 
the  tolerance  limits  of  5  km  and  5°  to  better  reflect  particularly  clean  or 
noisy  radar  data. 

A. 3  Spanning  over  Data  Gaps 

A  special  procedure  is  available  in  DRIVEA  to  selectively  override  obvious 
biases  and  extended  deviations  in  range,  azimuth  or  elevation  readings  by 
specifying  the  time  spans  over  which  it  is  desired  to  bridge  this  data  by 
smooth  spanning  cubic  polynomials.  This  function  is  implemented  following 
the  editor,  and  the  print-out  of  the  edited  data  includes  indication  of  these 
spanning  replacements  in  addition  to  the  data  added,  interpolated,  or 
modified.  Each  cubic  satisfies  readings  at  the  two  times  specified  plus  two 
readings  5  secs,  (and  no  less  than  2  secs.)  outside  the  span  of  discarded 
readings. 

The  user  must  specify  the  initial  (TI)  and  final  (TF)  times  over  which  one 
of  the  variables  (R,  A,  or  E)  is  to  be  spanned.  The  program  then  solves 
for  the  cubic  which  matches  values  of  this  variable  at  TI-5,  TI,  TF,  TF+5 
seconds.  Note  that  the  program  requires  data  to  be  available  for  at  least 
2  secs,  before  and  after  the  interval  to  be  spanned.  Otherwise  TI  is 
increased  and/or  TF  is  decreased  by  DRIVEA  to  satisfy  these  minimum 
requirements.  Finally  the  user  must  verify  that  the  variable  values  at  the 
four  specific  times  are  proper  for  generating  a  satisfactory  spanning  function. 
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Appendix  B 


Notation: 

P(k/n) 

P  ( k/N ) 

D  ( k ) 
ip  (k/n ) 

ip(k/N) 

<{»  (k+1  ,k) 


w(k) 

Q  ( k ) 

G(k) 

M(k) 


forward  filtering  position-velocity  vector  estimate  for  time  k 
based  on  observations  to  time  n. 

backward  filtering  final  estimate  for  position-velocity  vector, 
observed  data  vector  for  time  k. 

forward  filtering  error  covariance  matrix  associated  with  corres¬ 
ponding  estimate  for  P. 

backward  filtering  final  error  covariance  matrix  for  P(k/N). 

state  transition  matrix  from  K  to  k+1.  Thus,  when  purely 
integrating,  P{k+1)  =  4> (k+1  ,k >  *  P{k).  Note:  for  backward 
integration,  4»( k+1 )  =  (k+1  ,k)  is  obtained  by  using  -DT  in 

the  Runge-Kutta  calculations. 

covariance  matrix  of  the  estimated  radar  measurement  white  noise 
errors  at  time  k. 

covariance  matrix  of  the  unknown  factors  in  the  physical  model 
at  time  k. 

Kalman  filter  gain  for  forward  filtering  pass, 
weight  matrix  for  backward  filtering  pass. 


The  various  rotation  and  state  transition  matrices  are  necessary  to  connect 
the  trajectory  in  space  and  time,  but  the  noise  and  optimum  filtering 
estimation  is  described  below  ignoring  these  transformations.  The  operation 
is  also  best  understood  in  terms  of  a  single  rather  than  a  multi-variable 
process. 

Operation: 

During  the  forward  filtering  pass  the  error  covariance  matrices  <Kk+l,k) 
and  'Mk+l/k+l)  are  calculated  and  stored  for  time  k+1  from  time  k. 
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*(k+l/k)  =  *(k/k)  +  Q(k+1)  (B-l ) 

i.e.,  if  use  is  not  made  of  observations  to  update  the  estimated  position- 
velocity  vector  P,  this  estimate  will  be  poorer  with  time.  Note  that  except 
at  launch,  Q  is  a  few  orders  of  magnitude  less  than  <>. 

*  (k+l/k) 

G(k+1)  = -  (B-2) 

[♦(k+l/k)  +  W(k+1)] 

This  gain  matrix  determines  how  much  confidence  to  place  in  the  observed 
data.  If  the  radar  measurements  are  assumed  to  be  noisy  relative  to  the 
current  estimation  process,  reliance  on  each  data  point  should  be  reduced. 
If  physical  events  are  taking  place  which  are  unaccounted  for  in  the  model, 
or  if  the  radar  data  is  of  high  quality,  reliance  on  the  model  estimation 
should  be  reduced. 

P(k+l/k+l)  *  [1  -  G(k+1 )]  P(k)  +  G(k+1)  D(k+1)  (B-3) 

♦(k+l/k+1)  =  [1  -  G(k+1 )]  ♦(k+l/k)  (B-4) 

i.e.,  by  placing  some  reliance  on  the  recent  observation  the  estimated 
error  in  P  is  reduced. 

Besides  P(k+1),  i.e.  P(k+l/k+l),  both  ♦(k+l/k)  and  ♦(k+l/k+1)  are  stored 
on  a  file  for  later  use  on  the  backward  filtering  or  smoothing  pass. 

Since  forward  estimates  were  made  without  knowledge  of  future  observations, 
these  estimates  can  reasonably  be  expected  to  be  improved  by  a  "hindsight" 
backward  pass. 


♦(k/k) 

M(k)  =  (B-5) 

♦(k+l/k) 

This  weight  matrix  will  effectively  range  between  0  and  1  by  virtue  of 
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equation  (B-l). 


Let  y(k+l/k)  =  4>  -  P(k/k) 


Then 


P (k/N )  =  P(k/k)  +  M(k)[P(k+l/N)  -  y(k+l/k)]  (B-6) 

This  backward  filtering  adjustment  of  P  will  be  greater,  i.e.  more  influenced 
by  future  values,  when  M  approaches  1.  This  is  normally  the  case  in  mid¬ 
trajectory  when  the  error  covariance  matrix  is  not  changing  rapidly.  However, 
close  to  launch  the  state  vector  of  the .  rocket  is  known  reliably  and  the 
matrix  builds  up  from  zero.  Hence  on  the  backward  pass  M  approaches  0  near 
launch  and  the  trajectory  smoothly  assumes  the  original  forward  estimates 
P (k/k) . 

Finally, 


4»(k/N)  *  >Kk/k)  +  M(k)[  Kk+l/N)  -  ip(k+l/k)3  M(k)  T  (B-7) 

which  parallels  the  P(k/N)  logic  since  if  P (k/k)  is  greatly  influenced  by 
future  values  in  obtaining  P (k/N) ,  (k/N)  will  be  substantial ly  lower  than 

'/'(k/k).  Notice  that  'P ( k/N )  is  naturally  less  than  ip ( k/k)  and  is  equal 
only  for  k=N. 

In  practice,  the  error  covariance  matrix  that  is  computed  on  the  backward 
pass  sometimes  failed  to  remain  positive  definite.  This  problem  was  solved 
by  limiting  this  matrix  to  2%  of  the  forward  pass  error  covariance  matrix 
whenever  the  computed  value  for  a  diagonal  term  becomes  less  than  1%  of  the 
forward  estimate. 
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Appendix  C 


qualitative  Operation  of  the  Multi-Radar  Combining  Equations 
Notation: 

P(k)  Position-velocity  final  estimate  at  time  point  k. 

ip( k )  Final  error  covariance  matrix  associated  with  P(k) 

at  time  point  k 

Subscripts  1  and  2  will  identify  the  two  input  trajectory  estimates. 

The  composite  best  estimate  state  vector  is  given  by 

P(k)  =  [if'j-1  (k)  Px(k)  +  *2_1(k)  P2(k)]U1_1(k)  +  TP2"1(k)3"1  (C-l) 

The  associated  joint  error  covariance  matrix  is  given  by 

0(k)  =  [  i|'1"1(k)  +  *2”1(k)D"1  (C-2) 

Section  VI  of  Vol.  I  of  Reference  1  discusses  the  basis  of  the  original 
DRIVE4  multi-radar  estimating  program.  In  addition  to  the  primary 
equations,  C-l  and  C-2,  prior  weighting  and  scaling  was  available  to  de- 
emphasize  integration  results  further  away  in  time  from  filtered  radar 
data.  The  choice  of  weighting  was  subjective,  and  could  be  avoided  if  the 
error  covariance  matrices  reflected  the  confidence  in  the  estimate  of  the 
state  vectors  more  exactly.  Improvements  in  generating  error  covariance 
matrices  have  been  achieved  in  DRIVEB,  and  are  further  advanced  by  the  use 
of  backward  filtering.  The  desired  objective  of  using  the  above  equations 
exclusively  is  implemented,  with  the  exception  of  the  case  when  a  segment 
combines  radar  filtered  and  integrated  data?  in  this  case  the  radar 
filtered  results  are  copied  and  the  integrated  results  are  ignored. 
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Subroutine  COMBIN  is  called  during  each  time  step  to  generate  the  composite 
covariance  matrix,  state  and  acceleration  vectors  under  the  control  flag 
MODE.  The  reduced  applicable  MODE  values  are  described  below. 


M0D£=2  Trajectory  1  covariance  matrix,  state  and  accel¬ 
eration  vector  are  copied  to  obtain  the  composite 
estimate. 

M0DE=4  The  composite  trajectory  is  calculated  in  EVAL  using 
the  parameters  from  both  trajectories  according  to 
equations  C-l  and  C-2.  The  composite  acceleration 
estimate  is  obtained  by  using  the  same  relative 
weights  as  applied  to  position. 

M0DE=11  Trajectory  2  covariance  matrix,  state  and  accel¬ 
eration  vector  are  copied  to  obtain  the  composite 
matrix. 

Table  C  summarizes  the  mode  selection  algorithm  in  effect,  as  described 
above. 

The  condensed  procedure  described  above  has  not  been  extensively  tested, 
and  the  potential  exists  for  jumps  in  the  composite  state  vector  when 
switching  between  trajectories  1  and  2.  The  original  code  for  weighting 
and  scaling  in  COMBIN  using  other  MODE  values  has  therefore  been  retained, 
though  it  is  presently  non-functional.  As  noted  before,  the  desirable 
direction  is  to  have  the  separate  error  covariance  matrices  better  reflect 
the  results  of  filtering  and  integrating. 
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STATE  OF  TRAJECTORY  FILE  2 


Table  C 


;>Tni », 


TTiK 


*  Table 


STATE  OF  TRAJECTORY  FILE  1 


NULL 

RADAR 

INTEGRATING 

NULL 

NULL 

2 

2 

RADAR 

11 

4 

11 

INTEGRATING 

11 

2 

4 

Appendix  D 


r~ 


Modeling  of  the  Rocket  Dynamics 

The  modeling  of  the  rocket  dynamics  was  extended  considerably  from  the 
situation  when  only  after  burn-out  free  flight  trajectories  were  estimated 
in  DRIVE3.  A  straightforward  representation  was  developed  for  the  thrust, 
mass  loss,  and  drag  characteristics  of  most  of  the  rockets  encountered. 
Concurrently,  the  force  equations  were  developed  and  the  integration  procedure 
was  extended  to  include  the  thrusting  stages  from  take-off  to  burn-out. 

Assuming  thrusting  in  the  axial  direction  of  the  rocket  only,  the  thrust 
components  in  the  geocentric  inertial  x,  y,  z  system  are  given  by 

(o 

-  (C)  (B)  0 

where  the  axial  thrust  T  is  known,  B  rotates  this  thrust  to  the  local  CNR 
(east-north-vertical)  system  using  the  local  azimuth  Azr  and  elevation 
Elr  of  the  rocket  center-1  ine*,  C  rotates  the  ENR  thrust  components  of 
the  geocentric  x,  y,  z  values  by  the  right  ascension  «  and  declination  <5  of 
the  rocket  (see  Appendix  A,  Ref.  3). 

Adjusting  elevation  for  the  angle  of  attack  A  ,  the  azimuth  and  elevation 
of  the  rocket  is  determined  by  the  rocket  position  r  and  relative  velocity 
vector  vr. 

Azr  ,  El r  =  f(r,  Vp) 

vr  is  obtained  from  geocentric  quantities  and  the  earth  rotation  rate  w 
vr  =  7  -  "w  x  r 
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For  the  low  relative  velocities  near  lift-off,  the  launch  azimuth  and 
elevation  values  are  maintained  constant. 


D.l  Representation  of  the  Rocket  Characteristics 

A  piecewise  linear  approximation  is  made  for  each  significantly  different 
thrusting  phase.  Thrust  and  mass  are  expressed  as  straight  line  functions 
of  time  which  apply  until  the  next  time  phase  is  reached,  when  new  straight 
line  functions  apply.  For  the  drag  coefficient,  which  is  a  function  of 
Mach  number,  a  multiplier  is  applied  on  a  universal  Cq  curve  for  each 
time  segment.  This  universal  curve  is  given  by: 


CD  =  1,0  +  Vmach 


(for 


Vmach  — 


1) 


C 


D 


4.0 


mach 


+  1.0 


<f°r  vmach  >  '> 


For  a  particular  configuration  the  drag  coefficient  is  usually  10-25%  less 
during  the  power-on  stage  than  when  power  is  off. 

vmach  is  given  by  fVrl  /Vs  where  Vs  is  the  speed  of  sound.  For  the  U.S. 
standard  atmosphere,  1962,  the  speed  of  sound  is  given  as  a  function  of 
altitude  (Ref.  6).  This  relationship  is  reproduced  in  Figure  Dl. 


The  reference  cross-sectional  area  is  included  per  thrust  phase.  The 
thrusting  stage  (1,  2,  etc.)  is  also  given.  An  example  of  linearized 
characteristics  for  the  BLACK  BRANT  IV-B  (MOD  1)  without  payload  is  given 
in  Figure  D2.  This  rocket  is  represented  in  six  phases  as  follows: 


-46- 


Stage 

No. 


Time 

Secs. 

Thrust 

Lbs. 

AThrust 
per  sec. 

Weight 

Lbs. 

-  AWt. 

per  sec. 

Area 
Ft. 2 

CD 

Mult. 

Stag 

No. 

0 

27,500 

-333 

3294 

121 

0.56 

1.20 

1 

13.5 

23,000 

-5111 

1660 

35 

0.56 

1.20 

1 

18.0 

0 

0 

850 

0 

0.56 

.85 

2 

21.5 

8,400 

200 

850 

45 

0.56 

.70 

2 

31.0 

10,300 

-229 

420 

18 

0.56 

.70 

2 

35.5 

0 

0 

340 

0 

0.56 

.85 

2 

35.5  secs,  from  lift-off  the  rocket  is  in  free  flight. 

The  representation  of  the  thrust  history  needs  to  cover  pre-launch 
conditions  since  it  turns  out  in  many  instances  that  the  actual  thrusting 
and  lift-off  may  not  occur  until  a  few  seconds  after  launch  time.  This 
requirement  is  conveniently  implemented  by  providing  a  pre-launch  thrust 
equal  to  the  total  rocket  and  payload  weight,  i.e.  a  force  sufficient  to 
maintain  the  rocket  on  the  pad. 

D.2  Integration  Techniques 

Section  4.0  (Ref.  3)  details  the  adaptation  of  integration  techniques  to 
include  powered  flight.  Only  an  outline  of  the  approach  will  be  presented 
here. 


-49- 


The  Runge  Kutta  method  is  used  with  the  equations  of  motion  for  forward  or 
backward  time  integration. 


-  =  A  [p]  •  p  +  B  (p) 

dt 

where  "jr  is  the  position  vector 
TT  (p)  is  the  thrust  vector 
and  A  (p)  represents  the  free  flight  dynamics 

Then, 

Pn+1  =  Pn  +  l/6[Ki  +  2K2  +  2K3  +  £4] 

where 

AT  A[pn]  •  pn  +  aT  ’B(pn) 

I2  -  AT  A[pn  +  1/2  •  (pn  +  1/2  KX)  +  AT  B(pn  +  1/2  *1), 

K3  =  AT  A[pn  +  1/2  KZ]  *  (Pn  +  1/2  ^2)  +A  T  B(pn  +  1/2  K2), 

^4  «  AT  A[pn  +  I3]  •  (pn  +  K3)  +  AT  B(pn  +  K3), 

and  aT  is  the  sampling  interval  in  seconds.  Kj  is  evaluated  at 

time  tn,  K2  and  K3  at  tn  +  AT/2,  and  K4  at  tn  +  AT  or  tn+l* 

In  the  above  recurrence  relationships  if  the  thrust  matrix  IT  is  represented 
by  dividing  and  multiplying  by  the  position  vector  p,  the  resulting  normalized 
thrust  matrix  may  be  added  to  A(p).  A  procedure  identical  to  the  free-flight 
calculations  is  then  carried  out  to  give  the  desired  state  transition 


matrix  i>  such  that: 


Pn+i  =  (n+1  ,n)  pn 

Under  integration  or  filtering,  good  correspondence  to  expected  rocket 
performance  has  been  obtained  unless  the  rocket  malfunctioned  or 
experienced  large  wind  effects. 

D.3  Integration  Step  Size 

Some  raw  data  tapes  were  received  with  samples  at  1  second  rather  than 
the  usual  0.1  second  intervals.  It  was  found  that  if  the  integration  step 
size  was  allowed  to  increase  to  this  data  rate,  there  was  a  significant 
build-up  of  error  in  the  fourth  order  Runge-Kutta  integration  procedure. 
Step  size  less  than  0.1  second  resulted  in  negligible  change.  Subroutine 
RK  was  consequently  modified  to  loop  a  selectable  number  of  times  with  a 
specified  step  size.  The  maximum  step  size  is  set  at  0.1  second. 
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Appendix  E 


Least  Squares  Solution  for  an  Overdeternrined  Set  of  Equations 

The  classical  problem  in  trajectory  determination  when  ample  observations 
are  available  relative  to  the  variable  model  parameters  is  to  determine 
these  unknowns  so  as  to  minimize  the  sum  of  all  the  squares  of  the  separations 
between  the  observations  and  the  computed  trajectory  at  the  corresponding 
times.  These  unknowns  are  usually  the  initial  position  and  velocity  vectors 
plus  any  adjustable  model  parameters  such  as  drag  coefficient  and  geopotential 
terms  which,  once  fixed,  uniquely  determine  the  trajectory  to  within  the 
accuracy  of  the  model . 

This  optimization  problem  can  be  solved  directly  in  one  pass  if  the  process 
is  a  system  of  linear  equations.  The  theory  is  discussed  below.  However, 
the  process  is  generally  non-linear;  the  procedure  then  is  to  use  the  same 
theory  using  a  linearized  set  of  equations  for  the  partial s  of  the  variables 
(unknowns)  being  optimized.  Minimization  now  requires  an  iterative 
improvement  of  the  variables,  possibly  with  a  recomputation  of  the  partials  at 
the  observation  times  for  each  new  pass. 

E. 1  Linear  Equations 

Let  x  be  the  set  of  unknowns  to  be  estimated,  and  b  the  set  of  observations 
at  different  times.  At  each  sample  time  the  value  of  the  observation  may  be 
estimated  by  the  linear  equation  ax  where  a  is  a  fixed  set  of  coefficients 
corresponding  to  the  x's  and  the  observation  instant.  The  complete  set  of 
equations  is  thus  represented  by  Ax=b  where  each  line  of  this  equation  applies 
to  an  observation,  and  there  are  many  more  observations  (rows)  than  unknowns 
(columns). 

Each  observation  may  in  general  be  a  vector,  as  for  instance  in  trajectory 
determination  where  the  geocentric  rectangular  coordinates  are  usually 
taken  as  the  equally  weighted  metrics.  The  solution  follows  the  treatment 
of  Lanczos  (Ref.  5,  p.  156).  The  square  of  the  length  of  the  residual 
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vector  r  in  the  equation 


Ax  -  b  =  r  (E-l ) 

has  to  be  minimized. 

The  adjoint  matrix  A  is  obtained  by  transposing  the  rows  and  columns  of  A. 
For  example,  if 


-l  “ 

1  4 

A/ 

1  2  3 

A  = 

2  5 

then,  A  = 

3  6 

4  5  6 

_  *- 

The  solution  to  the  least  squares  problem  is  then  elegantly  given  by  the 
even  determined  system 

A  Ax  =  Ab  (E-2) 

which  has  just  as  many  equations  as  unknowns.  It  should  be  noted  that  the 
matrix  (A A)  is  symmetric  and  that  its  eigenvalues  are  positive  real. 

E.2  Non-linear  but  Differentiable  System 

Rather  than  express  results  in  the  general  notation  of  an  arbitrary  number 
of  unknowns,  the  procedure  is  outlined  below  for  three  unknowns  that  are  to 
be  optimized.  The  technique  will  be  seen  to  extend  without  restriction  to 
any  number  of  independent  unknowns. 

Assume  reasonable  original  estimates  xq  «(xi,  X2,  X3)  for  the  three 
unknowns.  Four  sets  of  equations  are  required,  one  for  the  system  and 
residuals  at  this  base  setting,  and  one  each  with  one  of  the  three 
variables  differentially  adjusted.  (This  is  the  simplest  form,  but  all 
three  variables  coula  be  varied  jointly  as  long  as  a  mutually  independent 
set  of  adjustments  is  achieved).  Let  these  adjustments  of  xj ,  X2,  X3 
one  at  a  time  be  given  by  axi,  ax 2,  AX3. 


The  residual  vector  for  n  observations  with  the  original  setting  is  ro  ,  and 
the  new  residuals  following  the  adjustments  are  given  by  rj,  r2,  r3  which 
may  be  determined  either  from  analytic  partials  of  the  system  model  or  from 
tests  on  the  model  with  the  adjustments  axj,  Ax2,  AX3  applied  one  by  one. 
Thus  we  have  the  vector 


*0 

xq  +  AX1 
X  =  XQ  +  A  X2 
x0  +  A  X3 

representing  four  trials  and  the  associated  residual  matrix 


01 

rll 

r21 

r31 

_ 

-  — 

- 

- 

- 

On 

rln 

r2n 

r3n 

for  n  observations  with  each  of  the  trials.  Expressed  in  q  units  of  x, 
where  we  wish  to  make  the  "observed"  residuals  (i.e.,  all  of  the  b's  of  Eqn. 
B-l )  zero,  the  problem  becomes  that  of  minimizing 


Rq  =  Aq 


(E-3) 


Hence  R  becomes  both  the  coefficients  of  the  unknowns  q  and  the  Ax-b  of 
equation  ( E - 1 ) -  The  procedure  for  minimizing  the  residuals  in  a  least 
squares  sense  is  then  suggested  by  equation  (E-2) 


R  R 


q0 


ql 

q2 

q-. 

w* 


=  0 


where  qj  is  in  units  of  and 


(E-4) 
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qQ  +  qj  +  q2  +  q3  =  1 


(E-5) 


R 


since  the  optimization  is  to  be  an  adjustment  about  Xq. 

The  four  equations  of  (E-4)  constitute  three  differential  relationships  such 
that  changes  to  xj,  X2,  X3  in  the  proper  combination  xj,  q2  4x3,  q3A  <3, 
will  minimize  the  residual  vector  ro.  The  solution  uses  the  property  that  the 
q^  are  in  the  ratio  of  the  row  sums  of  the  cofactors  of  R  R. 

In  practice  the  solution  may  be  implemented  by  obtaining  W  =  [R  R]~l,  the 
row  sums  wQ,  Wj,  w2,  w3>  and  their  total  T  =  wQ  +  Wj^  +  w2  +  w3  . 

Then  qj  =  Wj/T,  q2  =  w2/T,  q3  =  w3/T. 


An  iterative  procedure  may  be  incorporated  using  the  revised  values  for 
xg  as  the  base  point.  The  partials  determined  before  for  the  same 
observation  times  remain  applicable,  but  if  xq  underwent  a  large  change 
recomputation  of  the  partials  may  be  necessary.  It  should  be  noted  that 
monotonic  minimization  of  the  residual  metric  or  even  its  convergence  is 
not  assured  if  the  system  is  sufficiently  non-linear  and  if  the  original 
estimates  are  too  far  off  the  actual.  Criteria  for  stability  are  beyond 
the  scope  of  this  report. 

E.3  Example 

The  process  being  tested  is 

X  =  1.1A  (T  -  .98)0*85  +  o.gC 

Y  =  A/2.1  (T  -  .98)1*95  _  i.gc 

where  A,B,C  are  actually  2,1,1  and  X,Y  are  the  observations  at  any  time  T. 
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This  process  has  been  formulated  by  us  as 
X  =  A  (T-B)  +  C 

Y  =  A/2  (T-B)2  -  2C 

that  is,  the  expressional  form  is  approximately  known  and  the  associated 
optimum  values  for  A,B,C  are  to  be  determined  from  observations  on  the 
process. 

The  procedure  described  in  section  E.2  is  shown  in  operation  on  the  following 
pages.  A,B,C  are  originally  estimated  to  be  1,0,0  and  observations  are 
carried  out  at  times  T=l,2,3  and  4.  The  residual  vector  (Observed-Computed) 
for  X  and  Y  at  each  sample  time  is  thus  obtained.  The  residual  vectors  are 
also  obtained  for  the  cases  when  A,B,C  one  by  one  are  estimated  to  be  0.1 
larger.  This  gives  the  residual  R-matrix. 

The  square  root  of  the  sum  of  the  squares  of  the  residuals  before  each 
iteration  is  the  measure  of  the  closeness  of  the  estimate.  In  this  example 
the  original  estimate  and  sample  times  result  in  relatively  wide  variations 
of  the  residuals,  but  the  system  is  seen  to  approach  the  true  values  of 
A,B,C  after  three  iterations.  The  fourth  iteration  is  essentially 
superfluous.  The  optimum  obtained  for  A,B,C  is  1.78,  0.84,  1.02.  Further 
optimization  is  not  possible  unless  the  model  is  formulated  more  closely  to 
the  actual  process. 
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LEAST  SQUARES  OPTIMIZATION 
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LEAST  SQUARES  OPTIMIZATION 


ITERAT IO\ 
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LEAST  SQUARES  OPTIMIZATION  IT  r"?  AT  IOK  4 
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Appendix  F 


Formation  of  the  Driving  Noise  or  Uncertainty  Mat rix  Q 

Due  to  our  inexact  knowledge  of  the  physical  forces  influencing  the  rocke . 
trajectory,  the  state  vector  error  covariance  matrix  must  inevitably  increase 
with  time  in  the  absence  of  additional  radar  observations  for  updating  this 
vector.  This  phenomenon  is  mathematically  represented  by  adding  a  diagonal 
matrix  Q  to  the  error  covariance  matrix  at  each  step  of  the  integration 
process. 

In  the  original  DRIVE5  implementation  the  Q  element  corresponding  to  the  1th 
velocity  vector  was  given  by  Q(4-6)  =  (ACC(I)  •  QFACT  •  DT)2  where  ACC  is 
the  acceleration  at  the  time,  DT  is  the  integration  step  size,  and  QFACT  is  a 
judiciously  chosen  factor  usually  between  .01  and  .1.  The  corresponding 
position  vector  element  was  given  by  0(1-3)  =  Q (4-6)  •  DT2/4.  This  form 
makes  the  reasonable  assumption  that  model  uncertainties  are  greatest  at  times 
when  the  greatest  accelerations  are  experienced.  Operation  has  usually  been 
satisfactory  except  at  re-entry  when  model  and  data  sometimes  deviate. 

The  present  DRIVEB  in  the  normal  integration-filtering  mode  replaces  radar 
data  near  launch  by  an  optimum  trajectory  which  connects  smoothly  with 
subsequent  radar  data.  This  feature  plus  the  need  for  greater  emphasis  on 
data  at  re-entry  led  to  the  formulation  of  the  Ith  velocity  vector  as 
Q(4-6)  =  OVR  *[400(1)1*  (.01  •  QFACT  •  DT)2  where  0VR  is  the  vehicle  range 
at  the  time,  and  the  absolute  value  of  the  acceleration  is  used.  The  .01 
multiplier  gives  values  comparable  to  the  earlier  DRIVE5  version  for  20-30  km 
range  with  QFACT  relatively  unchanged.  The  position  vector  elements  are  as 
before  Q(l-3)  =  Q(4-6)  •  DT2/4.  This  form  has  proved  to  be  an  improvement. 
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Significance  of  Radar  Parameters  In  the  Filtering  Equations 

Radar  error  estimates  determine  the  measurement  error  covariance  matrix 
W(k )  which  in  turn  directly  affects  the  forward  pass  Kalman  filtering  gain. 
A  prime  concern  in  our  applications  is  to  obtain  the  best  transition  to 
radar  data  as  it  settles  down  from  launch  fluctuations  to  smoother  down 
range  readings. 

The  critical  choice  of  radar  parameters  for  most  rockets  is  therefore  where 
the  velocity  is  high  and  the  range  low,  and  where  the  data  is  first  considered 
reliable.  At  20  to  30  seconds  from  launch,  calculations  of  variances  for 
various  radars  have  established  the  following  factors  to  be  predominant 
over  other  noise  sources  by  two  or  more  orders  of  magnitude: 

Range  velocity  lag  coefficient  for  range  error 
Base  angular  variance  for  azimuth  and  elevation  error 
For  flights  where  the  radar  data  is  particularly  good  or  poor,  the  above 
parameters  may  be  appropriately  adjusted  from  nominal  to  better  represent 
the  operating  conditions  and  correspondingly  the  forward  filtering  gain. 
The  estimated  white  noise  matrix  W(k)  is  derived  from  the  variance  figures 
VARJRANGE,  AZ,  EL)  adjusted  for  range  and  elevation  of  observation,  and 
decomposed  into  the  geocentric  coordinate  system. 

Typical  range  error  variance  values  for  AN/FPQ-11  and  FPQ-6  radars  for 
a  rocket  range  rate  of  1.5  km/sec  at  20  to  30  secs,  from  launch,  and  a 
range  velocity  lag  coefficient  of  0.2  kyds/kyd/sec.  are  approximated  by 

(0.2  •  1.5)2  0.1  km2 

Corresponding  azimuth  and  elevation  error  values  for  base  angular  variance 
ranging  from  10.0  to  0.01  mil2  result  in  variance  figures  from  10”5  to  10"8. 
At  30  km  range  and  60°  elevation  these  are  equivalent  to  variances  of  0.005 
km2  and  below.  [10_5  .  302  .  cos{60°)]. 
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